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Resumo Estudamos problemas do calculo das varia(;6es e controlo optimo no 

contexto das escalas temporals. Especificamente, obtemos condi^oes 
necessarias de optimalidade do tipo de Euler-Lagrange tanto para la- 
grangianos dependendo de derivadas delta de ordem superior como para 
problemas Isoperimetricos. Desenvolvemos tambem alguns metodos dlrec- 
tos que permltem resolver determlnadas classes de problemas variacionais 
atraves de desigualdades em escalas temporals. No ultimo capftulo apresen- 
tamos operadores de diferen^a fraccionarios e propomos um novo calculo 
das variacoes fraccionario em tempo discreto. Obtemos as correspondentes 
condi^oes necessarias de Euler-Lagrange e Legendre, ilustrando depois a 
teoria com alguns exemplos. 

Palavras chave: calculo das varia0es, condi^es necessarias de op- 
timalidade do tipo de Euler-Lagrange e Legendre, controlo optimo, 
derivada delta, derivada diamond, derivada fraccionaria discreta, desigual- 
dades integrals, escalas temporals. 



Abstract We study problems of the calculus of variations and optimal control within 

the framework of time scales. Specifically, we obtain Euler-Lagrange type 
equations for both Lagrangians depending on higher order delta derivatives 
and isoperimetric problems. We also develop some direct methods to 
solve certain classes of variational problems via dynamic inequalities. In 
the last chapter we introduce fractional difference operators and propose 
a new discrete-time fractional calculus of variations. Corresponding 
Euler-Lagrange and Legendre necessary optimality conditions are derived 
and some illustrative examples provided. 

Keywords: calculus of variations, necessary optimality conditions of 
Euler-Lagrange and Legendre type, optimal control, delta-derivative, 
diamond-derivative, discrete fractional derivative, integral inequalities, time 
scales. 

2010 Mathematics Subject Classification: 26D10, 26D15, 26E70, 
34A40, 34N05, 39A12, 45J05, 49K05. 
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Introduction 



This work started when the author first met his advisor Delfim F. M. Torres in July, 2006. 
In a prehminary conversation Delfim F. M. Torres introduced to the author the Time Scales 
calculus. 

"A major task of mathematics today is to harmonize the continuous and the 
discrete, to include them in one comprehensive mathematics, and to eliminate 
obscurity from both." 



E. T. Bell wrote this in 1937 and Stefan Hilger accomplished it in his PhD thesis |60j in 
1988. Stefan Hilger defined a time scale as a nonempty closed subset of the real numbers 
and on such sets he defined a (delta) A-derivative and a A-integral. These delta derivative 
and integral are the classical derivative and Riemann integral when the time scale is the set 
of the real numbers, and become the forward difference and a sum when the time scale is 
the set of the integer numbers (cf . the definitions in Chapter [T|) . This observation highlights 
one of the purposes of Hilger's work: unification of the differential and difference calculus. 
But at the same time, since a time scale can be any other subset of the real numbers, an 
extension was also achieved (see Remark [321 pointing out this fact). As written in [23], the 
time scales calculus has a tremendous potential for applications. For example, it can model 
insect populations that are continuous while in season (and may follow a difference scheme 
with variable step-size), die out in (say) winter, while their eggs are incubating or dormant, 
and then hatch in a new season, giving rise to a nonoverlapping population. Motivated by 
its theoretical and practical usefulness many researchers have dedicated their time to the 
development of the time scales calculus in the last twenty years (cf. [Ml 125] and references 
therein). 

Since Delfim F. M. Torres works mainly in Calculus of Variations and Optimal Control 
and in 2006 only three papers |14^ \2T\ I61j . one of them by Martin Bohner [21], within this 
topic were available in the literature, we decided to invite Martin Bohner, whose knowledge 
in time scales calculus is widely known, to be also an author's advisor and we began to make 
an attempt to develop the theory. This was indeed the main motivation for this work. 

The two research papers [21] and [61] were published in 2004. In the first one, necessary 
conditions for weak local minima of the basic problem of the calculus of variations were es- 
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tablished, among them the Euler-Lagrange equation, the Legendre necessary condition, the 
strengthened Legendre condition, and the Jacobi condition, while in the second one, a cal- 
culus of variations problem with variable endpoints was considered in the space of piecewise 
rd-continuously A-differentiable functions. For this problem, the Euler-Lagrange equation, 
the transversality condition, and the accessory problem were derived as necessary conditions 
for weak local optimality. Moreover, assuming the coercivity of the second variation, a cor- 
responding second order sufficiency criterion was established. The paper [T3] was published 
in 2006 and there, were also obtained necessary optimality conditions for the basic problem 
of the calculus of variations but this time for the class of functions having continuous (nabla) 
V-derivative (cf . Definition E]) . The authors of [21 find this derivative most adequate to use 
in some problems arising in economics. 

In a first step towards the development of the theory we decided to study a variant of 
the basic problem of the calculus of variations, namely, the problem in which the Lagrangian 
depends on higher order A-derivatives. What seemed to be a natural consequence of the basic 
problem (as it is in the differential case) it turned out not to be! We only achieved partial 
results for this problem, being the time scales restricted to the ones with forward operator 
given by a{t) = ait + ao for some ai G M+ and ao € M, t S [a, p(6)]t (see Section 14.11 for 
the discussions). Another subject that we studied was the isoperimetric problem and here 
we were able to prove the corresponding Euler-Lagrange equation on a general time scale 
(cf. Theorem I63p . We were also interested in using direct methods to solve some calculus of 
variations and optimal control problems. This is particularly useful since even simple classes 
of problems lead to dynamic Euler-Lagrange equations for which methods to compute explicit 
solutions are very hard or not known. In order to accomplish this, some integral inequalitieqj 
needed to be proved on a general time scale. This is shown in Chapter [5l 

Another subject that interests us and that we have been working in is the development 
of the discrete fractional calculus [7B] and, in particular, of the discrete fractional calculus of 
variations. Recently, two research papers |15t [T6] appeared in the literature and therein it is 
suggested that in the future a general theory of a fractional time scales calculus can exist (see 
the book [77] for an introduction to fractional calculus and its applications). We have made 
already some progresses in this field (specifically for the time scale /iZ, h > 0), namely in the 
calculus of variations, and we present our results in Chapter [71 

This work is divided in two major parts. The first part has three chapters in which we 
provide some preliminaries about time scales calculus, calculus of variations, and integral 
inequalities. The second part is divided in four chapters and in each we present our original 
work. Specifically, in Section 14.11 the Euler-Lagrange equation for problems of the calculus 
of variations depending on the A-derivatives up to order r € N of a function is obtained. In 
Section [4.21 we prove a necessary optimality condition for isoperimetric problems. In Sections 



^We have also applied some integral inequalities to solve some initial value problems on time scales (cf. 
Chapter ini). 



5.11 and 15. 2| some classes of variational problems are solved directly with the help of some 
integral inequalities. In Section 16.11 some extensions of Gronwall's type inequalities are ob- 
tained while in Section [6. 2^ using topological methods, we prove that certain integrodynamic 
equations have C^[T] solutions. In Section [6.31 we prove Holder, Cauchy-Schwarz, Minkowski 
and Jensen's type inequalities in the more general setting of Oa-integrals. In Section F6.41 we 
obtain some Gronwall-Bellman-Bihari type inequalities for functions depending on two time 
scales variables. Chapter [7] is devoted to a preliminary study of the calculus of variations 
using discrete fractional derivatives. In Section [721 we obtain some results that in turn will 
be used in Section 17.31 to prove the main results, namely, an Euler-Lagrange type equation 
(cf. Theorem 1151 1) and a Legendre's necessary optimality condition (cf. Theorem I154|) . 

Finally, we write our conclusions in Chapter [8] as well as some projects for the future 
development of some subjects herein studied. 



Part I 



Synthesis 



Chapter 1 

Time Scales Calculus: definitions 
and basic results 



A nonempty closed subset of M is called a time scale and is denoted by T. The two most basic 
examples of time scales are the sets M and Z. It is clear that closed subsets of M may not be 
connected, e.g., [0, 1] U [2,3]. The following operators are massively used within the theory: 

Definition 1. The mapping a : T — > T, defined by a{t) = inf {s € T : s > t} with inf = 
supT (i.e., cr(M) = M if T has a maximum M) is called forward jump operator. Accordingly 
we define the backward jump operator p : T — )• T by p{t) = sup{s G T : s < t} with sup0 = 
inf T (i.e., p(m) = m liT has a minimum m). The symbol denotes the empty set. 

The following classification of points is used within the theory: A point t € T is called 
right-dense, right-scattered, left-dense and left- scattered if a{t) = t, a{t) > t, p{t) = t and 
p{t) < t, respectively. 

Let us define the sets T^ , inductively: T'' = T*^ = {t € T : t non-maximal or left-dense} 
and T'' =(¥" )'^, n > 2. Also, T^ = {t G T : t non-minimal or right-dense} and T^" 
= (T^u-i)«, n > 2. We denote T*^ n T« by T;J. 

The functions /x, i/ : T — )• [0, cxd) are defined by p{t) = a{t) — t and z^(t) = t — p{t). 



Example 2. If T = M, then a{t) = p{t) = t and p{t) = v{t) = 0. If T = [0, 1] U [2, 3], then 

a{t)- 



t if t G [0, 1)U[2,3]; 
2 ift = l, 



while, 



Also, we have that. 



P{t) 



Pit) 



t if t G [0, 1]U(2,3]; 

1 if t = 2. 

if t G [0, 1)U[2,3]; 

1 ift = l, 
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and, 

u(t) = l ° ifie[0,l]U(2,3]; 

For two points a, 6 € T, the time scales interval is defined by 

[a,b]j = {teT:a<t< b}. 

We now state the two definitions of differentiability on time scales. Throughout we will 
frequently write /'^(t) = /(cr(t)) and fP{t) = f{p{t)). 

Definition 3. We say that a function / : T — > M is A-differentiable at t G T*^ if there is a 
number f {t) such that for all e > there exists a neighborhood U oit (i.e., U = {t—5, t+5)r\T 
for some 5 > 0) such that 

\r{t) - /(s) - f''{t){<t) - s)\ < e\<y{t) - s\, for all s G U. 

We call f (t) the A-derivative of / at t. 

Remark 4. We note that if the number / (i) of Definition O exists then it is unique (see 
[581 [59], and also the recent paper [27| written in Portuguese). 

The A-derivative of order n G N of a function / is defined by / (i) = (/ " (t) ) , 
t G T''", provided the right-hand side of the equality exists, where / = /. 

Definition 5. We say that a function / : T ^ M is V -dijjerentiable at t G T^ if there is a 
number f {t) such that for all e > there exists a neighborhood U oit (i.e., U = {t—5, t+5)r\T 
for some (5 > 0) such that 

\fp{t) - f{s) - f^{t){pit) -s)\< e\pit) - s\, for ah s G U. 

We call /^(t) the V -derivative of / at t. 

Some basic properties will now be given for the A-derivative. We refer the reader to the 
works presented in [2^ for results on V-derivatives. 

Remark 6. Our results are mainly proved using the A-derivative. We point out that it is 
immediate to get the analogous ones for the V-derivative. This is done in an elegant way 
using the recent duality theory of C. Caputo 



Theorem 7. Jl^l Theorem 1.16] Assume / : T — > R is a function and let t G T^. Then we 
have the following: 

1. If f is A-differentiable at t, then f is continuous at t. 

2. If f is continuous at t and t is right-scattered, then f is A-differentiable at t with 



3. If t is right-dense, then f is A-differentiable at t if and only if the limit 

s~>-t S — t 

exists as a finite number. In this case, 

^ ^ s^t s-t 

4- If f is A-differentiable at t, then 

r{t) = f{t) + Kt)f^{t). (1.2) 

It is an immediate consequence of Theorem [7] that if T = M, then the A-derivative becomes 
the classical one, i.e., f^ = f while if T = Z, the A-derivative reduces to the forward 
difference f^{t) = Af{t) = f{t + 1) - fit). 



Theorem 8. |^ Theorem 1.20] Assume /, g : T ^ M are A-differentiable at t e T'^ . Then: 

1. The sum / + 5 : T ^ M is A-differentiable at t and (/ + g)^{t) = f^{t) + g^{t). 

2. For any number ^ G M, ^/ : T ^ M is A-differentiable at t and {if)^{t) = if^{t). 

3. The product fg:T^Ris A-differentiable at t and 

if9)^it) = fHt)9it) + r(i)5^(i) = fit)9^it) + fHt)g%t). (1.3) 

The next lemma will be used later in this text, specifically, in Section [4. II We give a proof 
here. 

Lemma 9. Let t € T*^ (t ¥" minTj satisfy the property p{t) = t < cr(i). Then, the jump 
operator a is not continuous at t and therefore not A-differentiable also. 

Proof. We begin to prove that lim^_j.(- a{s) = t. Let e > and take 5 = e. Then for all 
s & {t — 6,t) we have \cr{s) — t| < |s — t| < (5 = e. Since a{t) > t, this implies that a is not 
continuous at t, hence not A-differentiable by 1. of Theorem [71 D 

Now we turn to A-integration on time scales. 

Definition 10. A function / : T — ?> M is called rd-continuous if it is continuous at right-dense 
points and if the left-sided limit exists at left-dense points. 

Remark 11. We denote the set of all rd-continuous functions by Cj-d or CrdCT), and the set 
of all A-differentiable functions with rd-continuous derivative by C J^ or C J^ (T) . 
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Example 12. Consider T = Ufc^ol^^'^^ + ^]- -'^'-'^ ^^^^ time scale, 

' iftGUr=o[2A:,2A; + l); 



^^*^" I 1 iftGUr=o{2fc + l}- 



Joo 

lOO 



Let us consider t G [0, 1] PI T. Then, we have 



s^t t — s 



provided this hmit exists, and 



jA(l)^/(2)-/(l) 



2-1 ' 

provided / is continuous at t = 1. Let / be defined on T by 

r« if«€[0,l); 
I 2 ift> 1. 

We observe that at t = 1 / is rd-continuous (since hm(_^i- /(f) = 1) but not continuous (since 
/(I) = 2). 

Definition 13. A function F : T — )■ M is called an antiderivative of / : T'^ — > M provided 
F^{t) = f{t) for all f € T*^. 

Theorem 14. [24^, Theorem 1.74] Every rd-continuous function has an antiderivative. 

Let / : T'' :— )■ M be a rd-continuous function and let F : T — )• M be an antiderivative of /. 
Then, the A-integral is defined by jj f{t)At = F(r) — F(s) for all r, s € T. It satisfies 



rait) 

J /(r)AT = Mi)/(i), tGT\ (1.4) 



Theorem 15. f24^. Theorem 1.77] Let a, 6, c € T, ^ € M and f,g €Crd{T'^)- Then, 
rbt j./_,\ . /,M A^ rb j./,N A^ , rb 



1- /:[/(t) + <7(t)]Af = /:/(f)Af + /x^Ai. 

3. J^f{t)At = -£f{t)At. 

4. !'j{t)At = J:f{t)At + j'j{t)At. 

5. jy{t)At = 0. 

6. if\f{t)\<g{t) on[a,b]^, then 



f{t)At 

a 



10 



< / 9{t)At. 

J a 



One can easily prove j241 Theorem 1.79] that, when T = M then J^ f{t)At = f^ f(t)dt, 
being the right-hand side of the equahty the usual Rieniann integral, and when T = Z then 

jtfit)At = Etlfit)- 

Remark 16. It is an immediate consequence of the last item in Theorem 1151 that. 

f{t)<g{t), Vie[a,6]| implies [ f{t)At < f g{t)At. (1.5) 

Ja J a 

We now present the (often used in this text) integration by parts formulas for the A- 
integral: 

f r{t)g^{t)At = [(/5)(t)]*i - f f''{t)9{t)At- (1.6) 

J a J a 

' f{t)g^{t)At = [(/g)(t)]*i - t f^{t)g'^{t)At. (1.7) 



Remark 17. For analogous results on V-integrals the reader can consult, e.g., 

Some more definitions and results must be presented since they will be used in the sequel. 
We start defining the polynomials on time scales. Throughout the text we put No = {0, 1, . . .}. 

Definition 18. The functions gk,hk :T^ ^- M, defined recursively by 

go{t,s) = ho{t,s) = 1, s,teT, 

gk+i{t,s)= gk{cr{T),s)AT, hk+i{t,s)= /ifc(r,s)AT, A; G No, s,teT, 

J s J s 

are called the polynomials on time scales. 

Now we define the exponential function. First we need the concept of regressivity. 

Definition 19. A function p : T'' — ?> M is regressive provided 

1 + fi{t)p{t) + 

holds for all t G T'^. We denote by Tt the set of all regressive and rd-continuous functions. 
The set of all positively regressive functions is defined by 

7^+ = {p G 7^ : 1 + fi{t)p{t) > 0, for ah t G T''}. 

Tlieorem 20. 1251 Theorem 1.37] Suppose p G 7^ and fix to S T. Then the initial value 
problem 

y^=p{t)y, y{to) = l (1.8) 

has a unique solution on T. 

11 
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Definition 21. Let p G TZ and to S T- The exponential function on time scales is defined by 
the solution of the IVP (jl.Sp and is denoted by ep(-,to)- 

Remark 22. If T = M, the exponential function is given by 



with t, to £ IK and p : 



a continuous function. For T - 

t 
ep{t,to)= Y[[1+p{t)], 



the exponential function is 



T=to 

for t,tQ & 7j with to < t and p : Z — )■ M a sequence satisfying p{t) ^ — 1 for all t € Z. Further 
examples of exponential functions can be found in \2A\ Sect. 2.3]. 

The reader can find several properties of the exponential function in j24j . Let us now 
present some results about linear dynamic equations. 

For n G No and rd-continuous functions pj : T — ?> M, 1 < i < n, let us consider the nth 
order linear dynamic equation 



Ly = 0, where Ly = y^" + ^^pty^ 



(1- 



A function ?/ : T — )■ R is said to be a solution of equation (jl.Op on T provided y is n times 
delta differentiable on T" and satisfies Ly{t) = for all t G T'' . 

Lemma 23. \2^ p. 239] If z = {zi, . . . , z„) : T ^ M" satisfies for all t € T'' 

/ 1 ... \ 

: 1 ■•• : 

z'^ = A{t)z{t), where A= : •._•._ q (1.10) 

1 

\ -Pn -P2 -Pl J 

then y = zi is a solution of equation ()1.9p . Conversely, if y solves (|1.9p on T, then z = 
y,y'^, . . .,y'^"'') : T ^ M satisfies (fTTn]) for all t G T'"". 



Definition 24. [211 p. 239] We say that equation ()1.9p is regressive provided / + iJ.{t)A{t) is 
invertible for all t G T'', where A is the matrix in (jl.lOp . 

Definition 25. |24t p. 243] Let y/; : T — t- M be (m-1) times A-differentiable functions for 
all 1 < A; < 771. We then define the Wronski determinant W = W{yi, . . . ,ym) of the set 
{yi, • • • , Vm} by W{yi, ..., ym) = det V{yi, ..., ym), where 



/ yi 



V{yi,. 



vt 



y2 
y^ 



y-m 



\ 



\y^ 



m — 1 A rn — 1 






fl.ll) 
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Definition 26. A set of solutions {yi,...,y„,} of the regressive equation ()1.9p is called a 
fundamental system for ()1.9p if there is io S T'^ such that W{yi, . . . ,yn) 7^ 0. 

Definition 27. [Ml P- 250] We define the Cauchy function y : T x T'^" -> M for the linear 
dynamic equation ()1.9p to be, for each fixed s €T'^ , the solution of the initial value problem 

Ly = 0, y^'(a(s),s)=0, < i < n - 2, y'^"~\a{s),s) = I . (1.12) 

Theorem 28. 124[ p- 251] Suppose {yi, . . . ,yn} is a fundamental system of the regressive 
equation l[1.9\) . Let f G Crd- Then the solution of the initial value problem 

Ly = f{t), y^\to)=0, 0<i<n-l, 

is given by y(t) = f^ y{t,s)f{s)As, where y(t,s) is the Cauchy function for (|1.9p . 

We end this section enunciating three theorems, namely, the derivative under the integral 
sign, a chain rule and a mean value theorem on time scales. 



Theorem 29. j£^ Theorem 1.117] Let t^ e T'^ and assume A; : T x T ^ R zs continuous at 



{t,t), where i € T'^ with t > tg- -^'^ addition, assume that k{t, •) is rd-continuous on [tQ,a{t)]. 
Suppose that for each e > there exists a neighborhood U of t, independent of t ^ [tQ,a[t)], 
such that 

\k{a{t),T) - k{s,T) - k^^{t,T){a{t) - s)\ < e\a{t) - s\ for all s€U, 

where /c ^ denotes the delta derivative of k with respect to the first variable. Then, 

g(t) := / k(t, t) At implies g^{t) = / k^^{t,T)AT + k{a{t),t). 

Jto Jto 

Theorem 30. 124\ Theorem 1.90] Let / : M ^ M be continuously differentiable and suppose 
(7 : T — )• M is A- differentiable on T'^ . Then, f o g is A- differentiable and the formula 

if o g)^it) = 1^' f'igit) + hf,{t)g^it))dh^ g^it), t G T'^ 

holds. 

Theorem 31. 1251 Theorem I.I4] Let f be a continuous function on [a,b]f that is A- 
differentiable on [a, 6)t- Then there exist C,t € [a,b)j such that 

/A(o < m^ < /AM. 

b — a 
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Chapter 2 

Calculus of Variations and Optimal 
Control 



The calculus of variations deals with finding extrema and, in this sense, it can be considered a 
branch of optimization. The problems and techniques in this branch, however, differ markedly 
from those involving the extrema of functions of several variables owing to the nature of the 
domain on the quantity to be optimized. The calculus of variations is concerned with finding 
extrema for functionals, i.e., for mappings from a set of functions to the real numbers. The 
candidates in the competition for an extremum are thus functions as opposed to vectors 
in R", and this furnishes the subject a distinct character. The functionals are generally 
defined by definite integrals; the set of functions are often defined by boundary conditions 
and smoothness requirements, which arise in the formulation of the problem/model. Let 
us take a look at the classical (basic) problem of the calculus of variations: find a function 
y G C^la, b] such that 

C[y{-)] = [ Lit, yit), y'{t))dt -^ min, y{a) = ya, y{b) = yt, (2.1) 

Ja 

with a,b,ya,yb £ ^ and L{t,u,v) satisfying some smoothness properties. 

The enduring interest in the calculus of variations is in part due to its applications. We 
now present a historical example of this. 

Example 32 (Brachystochrones) . The history of the calculus of variations essentially begins 
with a problem posed by Johann Bernoulli (1696) as a challenge to the mathematical comunity 
and in particular to his brother Jacob. The problem is important in the history of the calculus 
of variations because the method developed by Johann's pupil, Euler, to solve this problem 
provided a sufficiently general framework to solve other variational problems [92]. 

The problem that Johann posed was to find the shape of a wire along which a bead initially 
at rest slides under gravity from one end to the other in minimal time. The endpoints of the 
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wire are specified and the motion of the bead is assumed frictionless. The curve corresponding 
to the shape of the wire is called a brachystochrone or a curve of fastest descent. 

The problem attracted the attention of various mathematicians throughout the time in- 
cluding Huygens, L'Hopital, Leibniz, Newton, Euler and Lagrange (see [92] and references 
cited therein for more historical details). 

To model Bernoulli's problem we use Cartesian coordinates with the positive y-axis ori- 
ented in the direction of the gravitational force. Let (a, ya) and (6, yb) denote the coordinates 
of the initial and final positions of the bead, respectively. Here, we require that a < b and 
ya < yh- The problem consists of determining, among the curves that have (a, ya) and (6, yh) 
as endpoints, the curve on which the bead slides down from (a, ya) to (6, yi,) in minimum time. 
The problem makes sense only for continuous curves. We make the additional simplifying 
(but reasonable) assumptions that the curve can be represented by a function y : [a,b] ^ M 
and that y is at least piecewise differentiable in the interval [a, b\. Now, the total time it takes 
the bead to slide down a curve is given by 

where / denotes the arclength of the curve, s is the arclength parameter, and v is the velocity 
of the bead s units down the curve from (a, b). 

We now derive an expression for the velocity in terms of the function y. We use the law 
of conservation of energy to achieve this. At any position {x,y{x)) on the curve, the sum of 
the potential and kinetic energies of the bead is a constant. Hence 

-mv'^{x) + mgy{x) = c, (2.3) 

where m is the mass of the bead, v is the velocity of the bead at (x, y(x)), and c is a constant. 
Solving equation ()2.3p for v gives 



m 



01= /'^^a... 

/2c 



Equality ()2.2p becomes 



'f-2<7y(x) 

We thus seek a function y such that T is minimum and y{a) = ya, y{b) = yb- 

It can be shown that the extrema for T is a portion of the curve called cycloid (cf. Example 
2.3.4 in Mi). 



Let us return to the problem given in ()2.ip and write the first and second order necessary 
optimality conditions for it, i.e., the well-known Euler-Lagrange equation and Legendre's 
necessary condition, respectively. 
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Theorem 33. (cf., e.g., ]92^ ) Suppose that L^u, L^v and L^^ exist and are continuous. Let 
y € C"'^[a,6] be a solution to the problem given in (|2.ip . Then, necessarily 

1. ■^L^{t,y{t),y'{t)) = Luit,y{t),y'{t)), te[a,b] (Euler-Lagrange equation) . 

2. L^^{t,y{t),y'{t)) > 0, t £ [a,b] (Legendre's condition). 

It can be constructed a discrete analogue of problem ()2.ip : Let a,b ^Tj. Find a function 
y defined on the set {a, a + 1, . . . , 6 — 1, 6} such that 

C[y{-)] = Y,L{t,y{t + l),Ky{t))^m\n, y{a) = y^, y{b) = yt, (2.4) 

t=a 

with ya,yb € M. For this problem we have the following result. 

Theorem 34. (cf., e.g., ^) Suppose that L^u, Luv and L^^ exist and are continuous. Let y 
be a solution to the problem given in (|2.4p . Then, necessarily 



1. ALi,{t) = Lu{t), i G {a, . . . , 6 — 2} (discrete Euler-Lagrange equation). 

2. Lyy{t) + 2Luv{t) + Luu{t) + Lyy{t + l) > 0, t G {o, . . . ,6 " 2}, (discrete Legendre's 
condition), 

where (t) = {t,y{t + 1), Ay{t)). 

In 2004, M. Bohner wrote a paper j21j in which it was introduced the calculus of variations 
on time scales. We now present some of the results obtained. 

Definition 35. A function / defined on [a, b]j x M is called continuous in the second variable, 
uniformly in the first variable, if for each e > there exists 6 > such that \xi — X2\ < S 
implies \f(t,xi) — f{t,X2)\ < e for all t E [a, 6]t- 

Lemma 36 (cf. Lemma 2.2 in [21]). Suppose that F{x) = J f{t,x)At is well defined. Lf fx 
is continuous in x, uniformly in t, then F'{x) = J fx{t,x)At. 

We now introduce the basic problem of the calculus of variations on time scales: Let 
a, 6 € T and L{t, u, v) : [a, b]j x M^ — )■ M be a function. Find a function y E Cjj such that 

rb 

C[y{-)] = / L{t,y''{t),y'^{t))At -^ min, y{a) = ya, y{b) = y^, (2.5) 

with ya,yb^ ^■ 

Definition 37. For / G C^.^ we define the norm 

11/11= sup ir(t)i+ sup i/^(t)|. 

A function y G Cj,(j with y{a) = ya and y{b) = yb is called a (weak) local minimum for problem 
()2.5p provided there exists 5 > such that C{y) < ^{y) for all y G Cj^j with y{a) = ya and 
y{b) = yb and \\y - y\\ < S. 
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Definition 38. A function t] € C^.^^ is called an admissible variation provided rj j^ and 
r]{a) = rj{h) = 0. 

Lemma 39 (cf. Lemma 3.4 in [2I])' Let y,ri ^ C^^ be arbitrary fixed functions. Put f{t,e) = 
L{t,y^{t) + erj'^{t),y^{t) + eri^{t)) and ^{e) = £{y + er]), e G M. If fe and fee are continuous 
in E, uniformly in t, then 

^\e)= [ [Lu{t,y^t),y^{t))v^t) + L,{t,y^t),y'^{t))rj^{t)]At, 



J a 

where [y](t) = (t,2/"(t),y^(t)). 

The next lemma is the time scales version of the classical Dubois-Reymond lemma. 
Lemma 40 (cf. Lemma 4.1 in [21]). Let g G Crd([a, 6]^). Then, 

b 

g{t)r]^it)At = 0, for all i] G Clij^{[a, 6]t) with ri{a) = r]{b) = 0, 

holds if and only if 

g{t) = c, on [a, b\j for some c ^M.. 

Nest theorem contains first and second order necessary optimality conditions for problem 
defined by ()2.5p . Its proof can be seen in |21j . 



Theorem 41. Suppose that L satisfies the assumption of Lemma\3{A If y & C^^ is a (weak) 
local minimum for problem given by ()2.5p . then necessarily 

1. L^[y\{t) = L„[y](t), t G [a, ^]x (time scales Euler-Lagrange equation). 

2. L^M{t) + l^{t){2LuM{t) + l^{t)Luu[y\{t) + (/«"W)*^™[y](^W)} > 0, t G [a,6]f' 
(time scales Legendre's condition), 

where [y](t) = {t,y''(t),y^{t)) and a* = ^ if a e M\{0} and 0* = 0. 

Example 42. Consider the following problem 

^[y{-)] = I [y'^{t)Y^t -^ min, y{a) = ya, y{b) = y^. (2.6) 



a 

2 



Its Euler-Lagrange equation is y {t) = for all t G [a, 6]| . It follows that y{t) = ct + d, and 
the constants c and d are obtained using the boundary conditions y{a) = ya and y{b) = yf,. 
Note that Legendre's condition is always satisfied in this problem since /z(t) > and L^u = 
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Chapter 3 



Inequalities 



Inequalities have proven to be one of the most important tools for the development of many 
branches of mathematics. Indeed, this importance seems to have increased considerably 
during the last century and the theory of inequalities may nowadays be regarded as an inde- 
pendent branch of mathematics. A particular feature that makes the study of this interesting 
topic so fascinating arises from the numerous fields of applications, such as fixed point theory 
and calculus of variations. 

The integral inequalities of various types have been widely studied in most subjects in- 
volving mathematical analysis. In recent years, the application of integral inequalities has 
greatly expanded and they are now used not only in mathematics but also in the areas of 
physics, technology and biological sciences [85]. Moreover, many physical problems arising in 
a wide variety of applications are governed by both ordinary and partial difference equations 
|64] and, since the integral inequalities with explicit estimates are so important in the study 
of properties (including the existence itself) of solutions of differential and integral equa- 
tions, their discrete analogues should also be useful in the study of properties of solutions of 
difference equations. 

An early significant integral inequality and certainly a keystone in the development of the 
theory of differential equations can be stated as follows. 

Theorem 43. If u ^ C[a,h] is a nonnegative function and 

u{t) < c+ du{s)ds, 

J a 

for all t G [a, h\ where c, d are nonnegative constants, then the function u has the estimate 

u{t) < cexp{d{t — a)), t G [a, b]. 

The above theorem was discovered by Gronwall [55] in 1919 and is now known as Gron- 
wall's inequality. The discrete version of Theorem 1431 seems to have appeared first in the work 
of Mikeladze [75] in 1935. 
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In 1956, Bihari [20] gave a nonlinear generalization of Gronwall's inequality by studying 
the inequality 

u{t)<c+ / dw{u{s))ds, 

J a 

where w satisfies some prescribed conditions, which is of fundamental importance in the study 
of nonlinear problems and is known as Bihari 's inequality. 

Other fundamental inequalities as, the Holder (in particular, Cauchy-Schwarz) inequality, 
the Minkowski inequality and the Jensen inequality caught the fancy of a number of research 
workers. Let us state and then provide a simple example of application (specifically to a 
calculus of variations problem) of the Jensen inequality. 

Theorem 44. [7H] If g G C{[a,b],{c,d)) and f G C{{c,d),M.) is convex, then 

J la9is)ds \ jtfigis))ds ^ 
\ b — a J b — a 

Example 45. Consider the following calculus of variations problem: Find the minimum of 
the functional C defined by 

>C[y(-)] = [\y'{t)]'dt, 
Jo 

with the boundary conditions y(0) = and 2/(1) = 1. Since f{x) = x^ is convex for all x € M 
we have that 



C[y{-)] > (^ y'{t)dt^ = 1, 



for all y G C^[0, 1]. Now, the function y{t) = t satisfies the boundary conditions and is such 
that C[y{t)] = 1. Therefore the functional C achieves its minimum at y. 

At the time of the beginning of this work there was already work done in the development 
of inequalities of the above mentioned type, i.e., Jensen, Gronwall, Bihari, etc., within the 
time scales setting (cf. [21 [831 EH [Ml [M])- We now state two of them being the others 
presented when required. The first one is a comparison theorem. 

Theorem 46 (Theorem 5.4 of [2]). Let a eT, y e Crd{T) and f € CrdiT"") and p G 7^+. 
Then, 

y'^it)<p{t)yit) + f{t), tGT^ 

implies 

y{t) < y{a)ep{t, a) + f ep{t, u{T))f{T)/\T, t G T. 

J a 

The next theorem presents Gronwall's inequality on time scales and can be found in [21 
Theorem 5.6]. 
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Theorem 47 (Gronwall's inequality on time scales). Let to € T. Suppose that u, a, b GCrrf(T) 
and b € 7^"^, 6 > 0. Then, 



implies 



u{t) < a{t) + / 6(r)u(r)AT for all t € T 

Jtn 



u{t)<a{t)+ a{T)b{T)eb{t,a{T))AT for all t G T. 
'to 



If we consider the time scale T = /iZ we get, from Theorem 1471 ^ discrete version of the 
Gronwall inequality (cf. [21 Example 5.1]). 

Corollary 48. // c, d are two nonnegative constants, a,b ^ hTL and u is a function defined 
on [a, b] n hlj, then the inequality 

u{t) <c+'^du{kh)h, t e[a,b]rihZ, 



h 



implies 



i-i 



t(t) < c+ ^cd{l + dh) ''h^^ h, t€[a,6]n/iZ. 



k=l 



Remark 49. We note that many new inequalities were accomplished by the fact that the proofs 
are done in a general time scale. For example, to the best of our knowledge, no Gronwall's 
inequality was known for the time scale T = IJfcgzt^'^ + 2]" ^'^^ ^^^® time scale, if a, 6 G T 
(a < b), then the A-integral is (see |58|): 



b M M-i 

fit)At= fit)dt+Yl 

''<' k=[a] 



^+h. 1 1 

f{t)dt + -f{k + -) 



+ / f{t)dt, 



where [t\ is the Gauss bracket. 
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Chapter 4 

Calculus of Variations on Time 
Scales 



In this chapter we will present some of our achievements made in the calculus of variations 
within the time scales setting. In Section I4.H we present the Euler-Lagrange equation for 
the problem of the calculus of variations depending on the A-derivative of order r € N of a 
function. In Section [4.21, we prove a necessary optimality condition for isoperimetric problems. 

4.1 Higher order A- derivatives 

Working with functions that have more than the first A-derivative on a general time scale 
can cause problems. Consider the simple example of f{t) = t^ . It is easy to see that / exists 
and f (t) = t + o"(t). However, if for example we take T = [0,1] U [2,3], then / doesn't 
exist at i = 1 because t + a{t) is not continuous at that point. 
Here, we consider time scales such that: 

(H) a{t) = ait + ao for some ai G M'^ and oq € M, t G [a, p{b)]j. 

Under hypothesis (H) we have, among others, the differential calculus (T = M, ai = 1, 
ao =0), the difference calculus (T = Z, ai = oq = 1) and the quantum calculus (T = {q : 
k G No}, with q > 1, ai = q, ao = 0). 

Remark 50. From assumption (H) it follows by Lemma [9] that it is not possible to have points 
which are simultaneously left-dense and right-scattered. Also points that are simultaneously 
left-scattered and right-dense do not occur, since a is strictly increasing. 

Lemma 51 (cf. [l5]). Under hypothesis (H), if f is a two times A-differentiable function, 
then the next formula holds: 

f''^{t) = aif'^''{t), teT^\ (4.1) 
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Proof. We have r^{t) = [f{t) + fi{t)f^{t)]'^ by formula ([L2D. By the hypothesis on a, 
function /i is A-differentiable, hence [f{t) + fi{t)f^{t)]^ = f^{t) + n^{t)f^''{t) + n{t)f^\t) 
and applying again formula ([OD we obtain f^{t) + fi^{t)f^''{t) + fi{t)f^\t) = /^'"(t) + 
fi^{t)f^''it) = (1 + /u^(t))/^'"(t). Now we only need to observe that n^{t) = a^{t) - 1 and 
the result follows. D 

We consider the Calculus of Variations problem in which the Lagrangian depends on 
A-derivatives up to order r G N. Our extension for problem given in ()2.5p is now enunciated: 



£[y(.)] = / L{t,y'^\t),y'^'- ^(t), . . . ,y-^''- (t),y^'-(t))At ^ min, 

J a 

y{a) = ya, y{p''~^{b))=yb, (P) 

y w = ya , y [p (^)) = yb > 

Assume that the Lagrangian L{t,uo,ui, . . . ,Ur) of problem ([P]) has (standard) partial 
derivatives with respect to uq, . . . ,Ur, r > 1, and partial A-derivative with respect to t of 
order r + 1. Let y G C^**, where 

C "^ = <y : [a,b]j ^ M. : y is continuous on T*^ > . 

Remark 52. We will assume from now on and until the end of this section that the time scale 
T has at least 2r + 1 points. Indeed, if the time scale has only 2r points, then it can be 
written as T = {a,a{a), . . . ,a^'-\a)}. Noting that p'-^b) = p'-^a^'"^ {a)) = a^'ia) and 
using formula (jl.ip we conclude that y{t) for t £ {a, . . . , a'''~^{a)} can be determined using the 
boundary conditions [y{a) = ya, ■ ■ ■ , y (a) = ya~^ and y{t) for t € {a^'{a), . . . , a'^^~^{a)} 
can be determined using the boundary conditions [y{p'''^^{b)) = yb,---,y {p'^'~^{b)) = 
y^~ ]. Therefore we would have nothing to minimize in problem (jPj) . 

Definition 53. We say that y^ € C^** is a weak local minimum for (jP|) provided there 
exists S > such that C{y^) < C{y) for all y € C^^ satisfying the constraints in ([P]) and 
\\y -y*||r-,oo < (5, where 

iiyll ._Vll7y«ll 
\\y\\r,oo ■ — / ^ iiy iicxDi 

j=0 
with yW = y'^"^" ' and ||y||oo := sup^g^fc'' |y(*)l- 

Definition 54. We say that rj € C^*" is an admissible variation for problem ([P]) if 

rj{a)=0, r^{p^-\b))=0 
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For simplicity of presentation, from now on we fix r = 3 (the reader can consult [71] for 
a presentation with an arbitrary r). 

Lemma 55 (cf. [IS]). Suppose that f is defined on [a,p^{b)]'f and is continuous. Then, 
under hypothesis (H), f^ f{'t)v^ {t)At = for every admissible variation rj if and only if 
f{t)=Oforallt€ [a,p^{b)]T. 

Proof. If f{t) = 0, then the result is obvious. 

Now suppose without loss of generality that there exists to ^ [a-, P^ (^)]t such that /(to) > 0. 
First we consider the case in which to is right-dense, hence left-dense or to = a (see Remark[SD]). 
If to = a, then by the continuity of / at to there exists a 5 > such that for all t € [tQ,to + 6)j 
we have /(t) > 0. Let us define i] by 

u) = / (* " *o)^^* " *o " ^^^ '^ * ^ [*0' *o + '^^TT' 
1 otherwise. 

Clearly rj is a C^ function and satisfy the requirements of an admissible variation. But with 
this definition for r] we get the contradiction 

rpHb) , /■to+<5 



/ f{t)ri^ (t)At = / f{tW {t)At > 0. 

J a J t{) 



'to 

Now, consider the case where to 7^ a. Again, the continuity of / ensures the existence of a 
6 > such that for all t € (to — 6, to -|- 6)j we have /(t) > 0. Defining r/ by 

ft) = / (* - *o + '^)^(* - *o - '^)^ if * e (to - <5,to + <5)t; 
1 otherwise, 

and noting that it satisfy the properties of an admissible variation, we obtain 

/ f{t)v'' (t)At = / f{t)v'' (t)At > 0, 

which is again a contradiction. 

Assume now that to is right-scattered. In view of Remark [501 all the points t such that 
t > to must be isolated. So, define rj such that rj'^ (to) = 1 and is zero elsewhere. It is easy 
to see that rj satisfies all the requirements of an admissible variation. Further, using formula 

dUD 

/ /(t)r?'^'(t)At = / /(t)77-'(t)At = A^(to)/(to)r?'^'(to) > 0, 

Ja Jt{) 

which is a contradiction. D 

Theorem 56 (cf. [12]). Let the Lagrangian L(t, uo, "Ui, ^2,^3) have standard partial deriva- 
tives with respect to uq,ui,U2,U3 and partial A-derivative with respect to t of order 4- On a 
time scale T satisfying (H), if y* is a weak local minimum for the problem of minimizing 

CP^iP) 

L{t,y-\t),y"'^{t),y"^\t),y'^\t)) At 
27 



CHAPTER 4. CALCULUS OF VARIATIONS ON TIME SCALES 



subject to 



y{a) = ya, y{p^(.b)) = y^, 

y^ia)=yly^{p\b))=yl 

y^\a) = yl y^'- {p\b)) = yl 



then y* satisfies the Euler-Lagrange equation 



L,, 



L 



«i< 



+ 1l 



U2 



Xl 



where (•) = {t,yfit),yf^{t),y-^\t),yfit)). 



3 u-i 



0, tG[a,/(6)]T, 



Proof. Suppose that y^, is a weak local minimum of C Let r] £ C^ he an admissible variation, 
i.e., r] is an arbitrary function such that r], tj and rj vanish at t = a and t = p'^{b). Define 
function $ : M ^> M by <I>(e) = £(y* + erf). This function has a minimum at e = 0, so we must 
have 

$'(0) = 0. 



(4.2) 



Differentiating $ with respect to e and setting e = 0, we obtain 

= j [Lu,{-)v-\t) + LuA.)r^-"\t) + L^,{.)r^-''\t)+L^^{-)^'^\t)]M. 
Since we will A-differentiate L„. , i = 1,2,3, we rewrite (|4.2p in the following form: 

= / {Lu,{-W\t) + LuA-W"''{t) + L^,{.)ri'^'^\t)+Lu,[-)v''\t)]/^t 

+ p{p\b)) [l^,^'^' + L^.rf'^ + L^.r?'^^^ + L^.r^"^' } {p\b)). (4.3) 
Integrating (|4.3|) by parts gives 



p3(b) 



{l„,(.)^'^' W - <(-)^'^' W - ^^.(O^'^^'^W - L^,{:)^^"''it)\ 



At 



+ 



LuA-W it) 



t=pHb) 



t=a 



+ [LuA-W''{t)tJa'''^+[LuA-)V^\t) 

+ p{p\b)) {L„„r?-' + L^,l^^''^ + Lu,r^^^' + L^,^"^'} {p\b)). 

2 

Now we show how to simplify ()4.4p . We start by evaluating rj'^ (a): 

77'^'(a)=r/'^(a)+M«)^-^(a) 

= ri{a) + p{a)7]^{a) + ^(a)air/^°'(a) 

= //(a)ai (^r?^(a) + //(a)?7^'(a)j 
= 0, 



t=p3(fe) 

t=a 



(4.4) 



(4.5) 
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where the last term of ()4.5p follows from ()4.ip . Now, we calculate rj'^ (a). By ()4.ip we have 
rj'^ {a) = aiT] (a) and applying ()1.2p we obtain 



air?^-(a) = ai (7?^(a) + ^(a)7?^'(a)) = 0. 



Now we turn to analyze what happens at t = p^{h). It is easy to see that if b is left-dense, then 
the last terms of (j4.4p vanish. So suppose that b is left-scattered. Since a is A-differentiable, 
by Lemma [9] we cannot have points that are simultaneously left-dense and right-scattered. 
Hence, p{b), p^{b) and p^{b) are right-scattered points. Now, by hypothesis, r] {p^{b)) = 0. 
Hence we obtain using (jl.ip that 

r,{p{b))-r^{pHb)) 
p{p^{b)) 

But Tj^p^ib)) = 0, hence r]{p(b)) = 0. Analogously, we have 

from what follows that r] {p[b)) = 0. This last equality implies 7/(6) = 0. Applying previous 
expressions to the last terms of (j4.4p . we obtain: 

r,-\p^{b))=v{pm=0, 

r^-\p^{b))=v{b)=0, 
^ ^^ ^ ^^ /^(p=^(&)) 

ri" ip(.b))-r (P' (b)) >?" (P' (&))->?" (P' (b)) 

P(p"(b)) " M(p"(fe)) 

/.(p3(6)) 

= 0. 

In view of our previous calculations, 



t=pHb) 
t=a 



+ p{p\b)) {K,r" + ^.x'?'^'^ + L,,,^""^' + L„3r?^^} (p3(6)) 

is reduced tqj 

L„3(p^(6))r?^'(p3(5)) + M(P'(6))L«3(P'W)^^'(P'W)- (4.6) 



^There is some abuse of notation: I/U3(p^(6)) denotes ^u3(-)lt=n3(6)' ''^^^ ^^' ^® substitute t in (■ 
{t,yf {t),yZ"^{t),yZ^\t),yf {t)) by p^(&). 
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Now note that 



V^'^p'm = v^\pHb))+f,{p'{b))r,^\p'{b)) 



and by hypothesis rj '^{p^{b)) = i] {p'^{b)) = 0. Therefore, 

l,{pHb))v^\p\b)) = -v^\p\b)), 
from which follows that (j4.6p must be zero. We have just simplified (j4.4p to 

= J^ {Luo{-K\t) - Lt^{-W\t) - L^^i-W'^^t) - Lt^{.)ii'^'%t)}^t. (4.7) 

In order to apply again the integration by parts formula, we must first make some transfor- 
mations in n" and rj °^. By (|4.ip we have 



ai 



and 

Hence, (j4.7p becomes 

rpHb) 



1 



aA^ 



r,^ -(i) = ^-- (t). 



(4i 



(4.9) 







L^^i.)rj'^\t) - L^^{.)rj'^\t) - ^ L^^{.)^^' ^ {t) - \L^^{.W'^\t)}M. (4.10) 



By the same reasoning as before, ()4.10p is equivalent to 

rpHb) 







Lu,{-)r^^\t) - L^^{-)r^"\t) - ^ L^^{.)^^' ^ {t) - \L^^{.)^^'^\t)]M 

(X\ Ci-\ ^ 



3 -A„<t3 1 rA a^A 1 rA„<TA2 



p{p\b))\L^,if -L^y--L^^v^ 



,2 "3 
^1 



and integrating by parts we obtain 

rp\b) 



= r {L^,{-)v''\t) - LtS-W\t) + ^L^l{.)^'^\t) + l<(-)^f^'^W}At 



1 



^^.(•)^'^ (i) 



t=P^(b) 



t=a 



1 



,2 "3 



^e3(-)^'^^W 



' A a» 1 






t=p4(b) 



t=a 
1 



(4.11) 



^^r^'^'^-^ii^'^^'K/^'W)- 



,2 «3 

^1 



Using analogous arguments to those above, we simplify (|4.1ip to 

fp\b) 



rp (b) , 11-1 

/ {L^,{-)rf\t) - L^^{.)rf\t) + ^L^^{-)rf' \t) + \Lf^{.)rf'^{t)Yt = 0. 

Ja ^ <1\ a-y ^ 
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Calculations as done before lead us to the final expression 



fp y") r q » Q 1 A 9 •? 1a-? •? -I 

/ {L^,{-)r (t) - L^S-K it) + -Lt i-W it) - -Lt i-)if (t)}At = 0, 



which is equivalent to 

/ L„„(.) - L,e,(-) + -<'(•) - -<(•) v''\t)At = 0. (4.12) 

Applying Lemma [551 to ()4.12p . we obtain the Euler-Lagrange equation 

^.o(-) - <(•) + ^L^^i-) - 3<(-) = 0' * e [a,p%b)]j. 

Oi\ Oil 

The proof is complete. D 

Following exactly the same steps of the proofs of Lemma [55] and Theorem [56] for an 
arbitrary r G N, one easily obtain the Euler-Lagrange equation for problem (|P]) . 

Theorem 57 (cf. |45j). (Necessary optimality condition for problems of the calculus of 
variations with higher-order A-derivatives) On a time scale T satisfying hypothesis (H), ify^ 
is a weak local minimum for problem (Q), then y^, satisfies the Euler-Lagrange equation 

^^(-irf-) ' <'(t,yf(t),yr'^(i),---,y*^^"W,yf(t))=o, (4.i3) 

i=o V«i/ ^ 

Remark 58. The factor ( ^ ) in (|4.13p comes from the fact that, after each time we apply 

the integration by parts formula, we commute successively a with A using ()4.ip [see formulas 
(|4.8p and (|4.9p ]. doing this X]!Ci J = 2 times for each of the parcels within the integral. 

Example 59. Let us consider the time scale T = {q^ : k € No}, a, 6 G T with a < b, and the 
problem 

^b(-)] = / iy^ it))' At -^ min, 

J a 

y(a) = 0, yipib)) = l, ^^"^^^ 

y^{a) = 0, y^{p{b)) = 0. 
We point out that for this time scale, a{t) = qt, p{t) = {q — l)t, and 

/^o = ^^('^l",(f ^ ljit)At= Y. ii-mit). 

t£[a,b)-j: 

By Theorem 1571 the Euler-Lagrange equation for problem given in (|4.14p is 

V(t) = o, te[a,p^{b)h. 
q 
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It follows that 

y^\t) =Ct + d, 

for some constants c, d € M. From this we get, using y^{a) = and y^{p{b)) = 0, 

y^{t) = c sAs-^ (t-a). 

J a p{o) - a 

Finally, using the boundary conditions y{a) = and y{p{h)) = 1, and defining 

["(') sAs 

Ja 



A 



B 



p{b) - a ' 



fr{i:rAr)As-Ajris-a)As 



P{b)i 



we get 



yit) = B 



I ( I tAt\ As-AJis- a)As 

Ja \J a / Ja 



t G [a, 6]t, 



(4.15) 



provided the denominator of B is nonzero. We end this example showing that this is indeed 
the case. We start by rewriting the denominator of B, which we denote by D, as 

r-p(b) 



D= i (r-ajAr) As 



Now, note that 



/-pW jf)(^s-a)As + a{pib)-a) f^W 
+ a / (s — a)As — r / (s — a)As. 

Ja pip) -a J a 



j U {T-a)AT\As = h3{p{b),a) 



Pib) 

{s - a)As = h2{p{b),a), 
where hi are given by Definition [181 It is known (cf. |241 Example 1.104]) that 



fc-i 



t-q^s 



h^{t,s) = YY 



, s,tGT. 



Hence, 



p_n PW-g'^Q I A p(fc) - g'-g n^=o g:ff + <P{^) - «) JL p(fc) _ g-^a 

From this it is not difficult to achieve 

^ ^ g(-/oW + q)(-p(^) + gQ)(« - P(^)g) 
(l + g)2(l + 5 + g2) 

Therefore, we conclude that D 7^ 0. We also get another representation for M.lSp . namely, 

(-t + a)(-t + ga)(a-tg) 



v{t) 



i-p{b) + a){-p{b) + qa){a - p{b)q) ' 



t G [a,b]T. 
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4.2 Isoperimetric problems 

In this section we make a study of isoperimetric problems on a general time scale (see [92] 
and [6] for continuos and discrete versions, respectively), proving the corresponding neces- 
sary optimality condition. These problems add some constraints (in integral form) to the 
basic problem. In the end, we show that certain eigenvalue problems can be recast as an 
isoperimetric problem. 

We start by giving a proof of a technical lemma. 

Lemma 60 (cf. [49j). Suppose that a continuous function / : T — ?• M is such that f'^{t) = 
for all t € T*^. Then, f{t) = for all t €T except possibly at t = a if a is right-scattered. 

Proof. First note that, since f'^{t) = 0, then /'^(t) is A-differentiable, hence continuous for all 
t € T''. Now, if t is right-dense, the result is obvious. Suppose that t is right-scattered. We will 
analyze two cases: (i) if t is left-scattered, then t ^ a and by hypothesis = f^ipit)) = f{t); 
(ii) if t is left-dense, then, by the continuity of /'^ and / at t, we can write 

Ve> 3(^1 >0: Vsi G (t - (^i, t]T, we have |r(si)-r(t)| <e, (4.16) 

Ve > 3(52 > : Vsa G (t - <52,i]T, we have |/(s2) - f{t)\ <e, (4.17) 

respectively. Let 5 = uiin{di,62} and take si € (t — 6,t)T. As cr(si) G (t — 5,t)T, take 
S2 = o-(si). By (|il6]) and ([iTT]) . we have: 

I - rit) + /(t)i = \risi) - nt) + /(t) - /(S2)i < \nsr) - f%t)\ + 1/(^2) - /(t)i < 2e. 

Since e is arbitrary, | — f'^{t) + f{t)\ = 0, which is equivalent to f{t) = f"{t). D 

We now define the isoperimetric problem on time scales. Let J : Cj^ — > M be a functional 
of the form 

J[yi-)]= f L(t,2/'^(t),y^(t))At, (4.18) 

J a 

where L(t, u, v) : [a, b]j x M x M — > M satisfies the assumption of Lemma [5U1 The isoperimetric 
problem consists of finding functions y G Cj,(j satisfying given boundary conditions 

y{a) = Va, y{b) = yb, (4.19) 

and a constraint of the form 

I[yi-)]= [ 9{t,y^t),y''{t))At = l, (4.20) 

where g{t, u, v) : [a, 6]^^ x R x M — >■ R satisfies the assumption of Lemma [39l and I is a specified 
real number, that takes ()4.18p to a minimum. 
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Definition 61. We say that a function y € Cj^ is admissible for the isoperimetric problem 
if it satisfies (fiA9]l and (020]). 

Definition 62. An admissible function y* is said to be an extremal for / if it satisfies the 
following equation: 



for all t G [a, b]j and some constant c. 

Theorem 63 (cf. |49j). Suppose that J has a local m,inim,um, at y^, € C^^ subject to the 
boundary conditions \4-l^ '^'^^ ^^c isoperimetric constraint Iji4.20\ ), and that y* is not an 
extremal for the functional I. Then, there exists a Lagrange multiplier constant A such that 
y-t satisfies the following equation: 

F,^{t,yUt),yHt))-Fu{t,yUt),yHt)) = 0, for all t e [a,b]^\ (4.21) 

where F = L — Xg. 
Proof. Let y* be a local minimum for J and consider neighboring functions of the form 

y = y* + eir]i + e2m, (4-22) 

where for each i € {1,2}, Si is a sufficiently small parameter (ei and £2 must be such that 
||y — ?/*|| < 6, for some 6 > 0), rii{x) G Cj^j and ?7j(a) = r]i{b) = 0. Here, rji is an arbitrary 
fixed function and rj2 is a fixed function that we will choose later. 

First we show that (j4.22p has a subset of admissible functions for the isoperimetric prob- 
lem. Consider the quantity 

rb 

nyi-)] = / g{t,yUt) + eiViit) + £2V2it),yHt) + eivHt) + £2V2it))^t. 

J a 

Then we can regard /[y(-)] as a function of ei and £2, say /[y(-)] = Q{£ii£2)- Since y* is 
a local minimum for J subject to the boundary conditions and the isoperimetric constraint, 
putting Q{ei,e2) = Q{£i,£2) — l we have that 

g(0,0) = 0. (4.23) 

By the conditions imposed on g, we have 



dQ '■'' 

deo 



(0,0) = / [guit,y^At),yHt))V2{t) + 9v{t,yUt),y^{t))r,iit)] At 

b 



gv{t,yUt),yHt))- I gu{r,yUr),yHr))^r 
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where (|4.24p follows from ()1.6p and the fact that r/2(a) = r]2{b) = 0. Now, consider the 
function 

E{t)=g,.{t,y^,{t),yt{t)) - [ 5«(T,yr(r), yf (r))Ar, t € [a,b]^. 

J a 

Then, we can apply Lemma HO] to show that there exists a function 772 S C^^ such that 



9vit,yUt),ytit)) - / gu{r,yUr),yHr))^r 



V^{t)At + 0, 



provided y* is not an extremal for /, which is indeed the case. We have just proved that 



dQ_ 



(0,0)/0. 



(4.25) 



Using (|4.23p and (j4.25p . the implicit function theorem asserts that there exist neighborhoods 
Ni and N2 of 0, Ni C {ei from (022])} n M and N2 C {eg from (022])} n M, and a function 
£2 : Ni ^fM. such that for all £1 G A'^i we have 



Q{£l,£2{£l)) = 0, 

which is equivalent to Q{£i,£2{£i)) = ^- Now we derive the necessary condition (j4.2ip . 
Consider the quantity J{y) = K{£i,£2). By hypothesis, K is minimum at (0,0) subject to 
the constraint Q(0,0) = 0, and we have proved that VQ(0,0) 7^ 0. We can appeal to the 
Lagrange multiplier rule (see, e.g., [92^ Theorem 4.1.1]) to assert that there exists a number 
A such that 

V(ir(0,0)-AQ(0,0)) =0. (4.26) 

Having in mind that 7?i(a) = r/i(6) = 0, we can write: 



||(0, 0) = y [L„(t, yl{t),yt{t)Wi{t) + L,(t, y:{t),yt{t))rjt{t)\ Ai 



J a 



L,{t,y:it),y^it))- / L„(r,y,^(T),yf(r))AT 



r?f(i)At. 



(4.27) 



Similarly, we have that 

rb 



S'->^r 



9vit,yUt),y^{t))- / g4T,yUr),ytiT))AT 



r]i{t)At. (4.28) 



Combining (g^S]), (ji:^ and (|05|) . we obtain 



L„(-) 



L„(-)Ar-A 



5.(-)-y gJ")ArHr?f(t)At = 0, 



where (■) = {t,y^{t),y^{t)) and (••) = (T,<(T),yf (r)). Since r?i is arbitrary. Lemma HO] 
implies that there exists a constant d such that 

Lv{-) - Xgvi-) - (J [Lu{-) - A5„(")]Ar j = d, tG [a, 6]^, 



35 



CHAPTER 4. CALCULUS OF VARIATIONS ON TIME SCALES 



or 



F^,{-) - f F,,{-)At = d, (4.29) 

J a 

with F = L — Xg. Since the integral and the constant in (|4.29p are A-differentiable, we obtain 
the desired necessary optimahty condition (|4.2ip . D 

Remark 64. Theorem [63] remains vahd when y^ is assumed to be a local maximizer of the 
isoperimetric problem ()4.18p - ()4.20p . 

Remark 65. In Theorem 1631 we assume that y* is not an extremal for /. This corresponds to 
the normal case. A version of Theorem 1631 that involves abnormal extremals can be found in 
|1U] . A different hypothesis excluding abnormality is given in 



Example 66. Let a, —a € T and suppose that we want to find functions defined on [—a, a\j 
that take 

J[yi-)] = r y'-itW 

J —a 

to its largest value (see Remark I64p and that satisfy the conditions 

y{-a) = via) = 0, I[y{-)] = f ^l + (y^{t)YM = I > 2a. 

Note that if y is an extremal for /, then y is a line segment |21j . and therefore y{t) = 
for all t E [— a, ajf. This implies that I(y) = 2a > 2a, which is a contradiction. Hence, I 
has no extremals satisfying the boundary conditions and the isoperimetric constraint. Using 
Theorem ESI let 



F = L - A5 = y" - A^l + {y^y . 
Because 



y^ 

Fx — 1, F^ — A- 



Vi + W? 

a necessary optimality condition is given by the following dynamic equation: 

X\ , ^ =1 -1 = 0, t(^[-a,a]^\ 

If we restrict ourselves to times scales T satisfying hypothesis (H) of Section [4. II it follows 
that the same proof as in Theorem [63] can be used, mutatis m.utandis, to obtain a necessary 
optimality condition for the higher-order isoperimetric problem (i.e., when L and g contain 
higher order A-derivatives) . In this case, the necessary optimality condition ()4.2ip is gener- 
alized to 



B-..1il 


2 

F; 


again F = L — Xg. 





{t,yf{t),yf"^{t),...,y^,^''-\t),yf{t))=0, 
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4.2.1 Sturm— Liouville eigenvalue problems 

Eigenvalue problems on time scales have been studied in [3]. Consider the following Sturm- 
Liouville eigenvalue problem: find nontrivial solutions to the A-dynamic equation 

y^'(i) + <?(%"(*) + Ay" (t) = 0, te[a,b]^\ (4.30) 

for the unknown y : [a, b]j —^ M subject to the boundary conditions 

y{a) = y{b) = 0. (4.31) 

Here q : [a, 6]| — >■ R is a continuous function. 

Generically, the only solution to equation ()4.30p that satisfies the boundary conditions 
()4.3ip is the trivial solution, y{t) = for all t E [a,6]Tr- There are, however, certain values 
of A that lead to nontrivial solutions. These are called eigenvalues and the corresponding 
nontrivial solutions are called eigenf unctions. These eigenvalues may be arranged as — oo < 
Ai < A2 < . . . (cf. [31 Theorem 1]) and Ai is called the first eigenvalue. 

Consider the functional defined by 

J[y{-)] = l\{y''?{t) - q{tWf{t)W, (4.32) 



and suppose that y^ G C^^ (functions that are twice A-differentiable with rd-continuous second 
A-derivative) is a local minimum for J subject to the boundary conditions (|4.31|) and the 
isoperimetric constraint 

m-)]= f\yn\t)At = i. (4.33) 

Ja 

If y* is an extremal for /, then we would have —2y"[t) = 0, t S [a, h\j. Noting that y[a) = 0, 
using Lemma [60] we would conclude that y{t) = for all t G [a, ^Jt- No extremals for / can 
therefore satisfy the isoperimetric condition (|4.33p . Hence, by Theorem [63l there is a constant 
A such that y* satisfies 

F^{t,y'i{t),y^{t)) - F^{t,y'i{t),y^{t)) = 0, t € [a,6]|', (4.34) 

with F = v"^ — qu^ — \u^. It is easily seen that ()4.34p is equivalent to ()4.30p . The isoperimetric 
problem thus corresponds to the Sturm-Liouville problem augmented by the normalizing 
condition (|4.33p . which simply scales the eigenfunctions. Here, the Lagrange multiplier plays 
the role of the eigenvalue. 

The first eigenvalue has the notable property that the corresponding eigenfunction pro- 
duces the minimum value for the functional J . 



Theorem 67. Let Ai he the first eigenvalue for the Sturm-Liouville problem { 4-^0\ ) with 



boundary conditions ^.31 ), and let yi be the corresponding eigenfunction normalized to satisfy 



the isoperimetric constraint ^.33 ). Then, among functions in C^^ that satisfy the boundary 



conditions 1^4-31 ) and the isoperimetric condition ^.33^ , the functional J defined by ^ 
has a minimum at yi. Moreover, J[yi{-)] = Ai. 
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Proof. Suppose that J has a minimum at y satisfying conditions (|4.3ip and (|4.33p . Then y 
satisfies equation ()4.30p and multiplying this equation by y'^ and A-integrating from a to p{b) 
(note that y is defined only on [a, 6]| ), we obtain 

rPib) .. rp{b) rp{b) 

/ y-(t)y^ (t)At + / q{t){y-f{t)^t + A / (y'^)2(t)At = 0. (4.35) 

Ja Ja Ja 

Since y{a) = y{b) = 0, we get 

fpib) ,,> rp(b) 

/ y'^{t)y^\t)M = [y{t)y\t)\tT - / {y'')\t)M 

J a J a 

= yip{b))y^ipib))- / (2/^)'(t)At. (4.36) 

J a 

If b is left-dense it is clear that (|4.36p is equal to — J (y )^(t)At. If b is left-scattered (hence 
p{b) is right-scattered) we obtain for ()4.36p 



yipmy'^iPib)) - I (y^)'(t)At = ^(p(^)) ^.(^(^j) - / (y'')!*)^^ 



a 



.(P(^))^ r^'(,A)^(OA. 



\y^)\t)At, 



where the last equality follows by p.4p . With the help of (|1.4p it is easy to prove that 
jf^qmynHt)At = jlq{t){y^)\t)M and j^^"^ {yn\t)At = j'^{ynHt)At. By (gSSl), 
()4.35p reduces to 

\{y^f{t)-q{t){yn\t)]At = \, 

that is, J[2/(-)] = ^- Due to the isoperimetric condition, y must be a nontrivial solution to 
(j4.30p and therefore A must be an eigenvalue. Since there exists a least element within the 
eigenvalues, Ai, and a corresponding eigenfunction yi normalized to meet the isoperimetric 
condition, the minimum value for J is Ai and J[yi(-)] = Ai. D 

4.3 State of the Art 

The results of this chapter are already published in the following international journals and 
chapters in books: \i5\ I49j . In particular, the original results of the papers [451 H9] were 
presented in the International Workshop on Mathematical Control Theory and Finance, Lis- 
bon, 10-14 April, 2007, and in the Conference on Nonlinear Analysis and Optimization, June 
18-24, 2008, Technion, Haifa, Israel, respectively. It is worth mentioning that nowadays other 
researchers are dedicating their time to the development of the theory of the calculus of 
variations on time scales (see [lOl [231 E21 EH UM EH [96] and references therein). 
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Chapter 5 

Inequalities applied to some 
Calculus of Variations and Optimal 
Control problems 



In this chapter we prove a result that complements Jensen's inequality on time scales (cf. 
Proposition [5U|) and state some useful consequences of it. These are then applied in Section [^T^ 
to solve some classes of variational problems on time scales. A simple illustrative example is 
given in Section 15.2.11 

The method here proposed is direct, in the sense that permits to find directly the optimal 
solution instead of using variational arguments and go through the usual procedure of solving 
the associated A Euler-Lagrange equation. This is particularly useful since even simple 
classes of problems of the calculus of variations on time scales lead to dynamic Euler-Lagrange 
equations for which methods to compute explicit solutions are not known. A second advantage 
of the method here promoted is that it provides directly an optimal solution, while the 
variational method on time scales is based on necessary optimality conditions, being necessary 
further analysis in order to conclude if the candidate is a local minimizer, a local maximizer, 
or just a saddle. Finally, while all the previous methods of the calculus of variations on time 
scales only establish local optimality, here we provide global solutions. 

The use of inequalities to solve certain classes of optimal control problems is an old idea 
with a rich history [30l [38l [56l EH [89] . We trust that the present study will be the beginning 
of a class of direct methods for optimal control problems on time scales, to be investigated 
with the help of dynamic inequalities. 

5.1 Some integral inequalities 

The next theorem is a generalization of the Jensen inequality on time scales. 
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Theorem 68 (Generalized Jensen's inequality [94j). Let a,6 G T and c,d S R. Suppose 

f : [a, b]j — 7> (c, d) is rd-continuous and F : {c, d) ^ M is convex. Moreover, let h : [a, b]j -^ M 
be rd-continuous with 



\h{t)\At >0. 

J a 

Then, 



j^\hit)\Fim)At ^^ nt\hit)\mAt \ 

Ia\Ht)\^t [ fjh{t)\At I' 



The next observation is crucial to solve variational problems. Follows the statement and 
a proof. 



Proposition 69 (cf. j22j). If in Theorem [23 F is strictly convex and h{t) ^ for all 
t € [a,^]^^, then the equality in i5.1\) holds if and only if f is constant. 

Proof. Consider xq G (c, d) defined by 

J^\h{t)\f{t)At 

From the definition of strictly convexity, there exists in gM. such that 

F{x) — F{xq) > m{x — xq), 

for all X G (c, (i)\{xo}. Assume / is not constant. Then, /(to) 7^ ^0 for some to € [a,^]^^. We 
split the proof in two cases, (i) Assume that to is right-dense. Then, since / is rd-continuous, 
we have that /(t) ^ xo on [to, to + 5)t for some 5 > 0. Hence, 



\h{t)\F{f{t))At- f \h{t)\AtFixo) = I \hit)\[Fif{t))-Fixo)]At 



>m f |/i(t)|[/(t)-xo]At 

J a 



0. 
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(ii) Assume now that to is right-scattered. Then [note that f^ ° f{t)At = ^J.{to)f{to)], 
'\h{t)\F{f{t))At- f \h{t)\AtF{xo) 

J a 

\hmF{m) - F{xo)]At 

'a 

'"\hmF{f{t))-F{xo)]At+ r'" \hmF{f{t)) - F{xo)]At 

J to 

+ f \hmFum-F{x,)w 

J(T(to) 

> j'\h{t)\[F{f{t))-F{xQ)]At + m r ° \h{t)\[f{t)-xo]At 
+ f \h{t)\[F{f{t))-F{x,)]At 

Ja(to) 

>m\ r\hmf{t)-xo]At+ r*" \hmf{t)-xo]At 

\ J a J to 



+ I Ht)\[f{t)-x,]At\ 

Ja(to) ) 

= m f \hmf{t)-xo]At = 0. 

J a 

Finally, if / is constant, it is obvious that the equality in ()5.ip holds. D 

Remark 70. If in Theorem 1681 F is a concave function, then the inequality sign in ()5.ip must 
be reversed. Obviously, Proposition [69] remains true if we let F to be strictly concave. 

Before proceeding, we state a particular case of Theorem 



Corollary 71. Let a = q"' and b = q'^ for some n,r?i € No with n < m. Define f and h on 
j^n^^m-ij ^y^^ assume F : (c, d) — )■ M is convex, where {c,d) D [/(<?")> /('?"'~"^)]„no . // 

rra— 1 



J2(^-l)q'\h{q')\>0, 



k=n 

then 



j:T=nq'\Hq')\F{f{q'')) ^ ^ ( Ek=n q'\Kq')\f{q') 



ET=nq'\Hq')\ ~ V T.k=nq'Hq')\ J' 

Proof. Choose T = q^o = {^fc . j. ^ f^^}^ ^ > i^ i^ Theorem EH □ 



We now present Jensen's inequality on time scales, complemented by Proposition 169 
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Theorem 72 (cf. |22j). Let o, 6 G T and c, d € M. Suppose f : [a, b]j — ?■ (c, d) is rd-continuous 
and F : (c, (i) — )> M is convex (resp., concave). Then, 

(resp., the reverse inequality) . Moreover, if F is strictly convex or strictly concave, then 
equality in i5.2\) holds if and only if f is constant. 

Proof. Particular case of Theorem 1681 and Proposition [SU] with h{t) = 1 for ah t € [a,b]j. D 

We now state and prove some consequences of Theorem [72l 

Corollary 73 (cf. [22]). Let a,b e T and c,d G M. Suppose f : [a, 6]f. -^ {c,d) is rd- 
continuous and F : (c, d) — >■ M is such that F" > (resp., F" < 0). Then, 

(resp., the reverse inequality). Furthermore, if F" > or F" < 0, equality in i5.3\) holds if 
and only if f is constant. 

Proof. This follows immediately from Theorem [72] and the facts that a function F with F" > 
(resp., F" < 0) is convex (resp., concave) and with F" > (resp., F" < 0) is strictly convex 
(resp., strictly concave). D 

Corollary 74 (cf. [22]). Let a,b e T and c,d G M. Suppose f : [a, 6]f. -)■ {c,d) is rd- 
continuous and ip,ip : {c,d) — )■ M are continuous functions such that ip^^ exists, tp is strictly 
increasing, and ip o ip~^ is convex (resp., concave) on Ini(ip). Then, 

^-1 / LNX/(t))AA ^ ^_, / />(/(t))At \ 

(resp., the reverse inequality). Furthermore, if ip ° V~^ ^-^ strictly convex or strictly concave, 
the equality holds if and only if f is constant. 

Proof. Since if is continuous and (p o f is rd-continuous, it follows from Theorem [72] with 

f = ip o f and F = ip o ip~^ that 

J>o^-l)((v.o/)(t))At ^ ^ / />o/)(OAA 

b-a -^^ ^ ' y b-a J ' 

Since ip is strictly increasing, we obtain 

^-1 / />(/(t))At \ ^ / />(/(t))At \ _ 

I 6 — a i I 6 — a i 

Finally, when -0 o (^~^ is strictly convex, the equality holds if and only if 99 o / is constant, or 
equivalently (since (p is invertible), / is constant. The case when ipoip~^ is concave is treated 
analogously. D 
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Corollary 75 (cf. [22]). Assume f : [a,b]j —^Mis rd- continuous and positive. If a < or 
a > 1, then 

f\f{t)rM>{b-af~'^(ff{t)M 

J a \J a J 

IfO<a<l, then 

f\f{t)rAt<{b-a)'~''( f'fit)At 

Ja \Ja J 

Furthermore, in both cases equality holds if and only if f is constant. 
Proof. Define F{x) = x", x > 0. Then 

F"(x) = a(a - l)x"-^ x > 0. 



Hence, when a < or a > 1, F" > 0, i.e., F is strictly convex. When < q < 1, F" < 0, 
i.e., F is strictly concave. Applying Corollarv 1731 with this function F, we obtain the above 
inequalities with equality if and only if / is constant. D 

Corollary 76 (cf. [22]). Assume f : [0,^*]^ -^M is rd-continuous and positive. If a < —1 or 
Q > 0, then 



fwA/jf'"''^'"^''-'- 



l+a 



If -1 < a <0, then 

Furthermore, in both cases the equality holds if and only if f is constant. 

Proof. This follows from Corollary [7S] by replacing / by 1// and a by —a. D 

Corollary 77 (cf. [2^). // / : [a, &]| — )• M is rd-continuous, then 

f e^(*)At >{b- a)e^^-'^(*)'^* . (5.4) 

Ja 

Moreover, the equality in ( [5.^[ ) holds if and only if f is constant. 

Proof. Choose F{x) = e^, x e M, in Corollarv 1731 D 

Corollary 78 (cf. [22]). If f ■ [o-yb]!^ — >■ M is rd-continuous and positive, then 

^ ln(/(t))At <{b-a) In (^^ ^ f{t)At^ . (5.5) 

Moreover, the equality in i5.5\) holds if and only if f is constant. 

Proof. Let F{x) = ln(x), x > 0, in Corollarv 1731 D 
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Corollary 79 (cf. [22]). // / : [a, bj^f —^R is rd-continuous and positive, then 

^ /(t)ln(/(t))At>^ /(t)Atln(^^^ /(t)At) . (5.6) 

Moreover, the equality in (j5.6p holds if and only if f is constant. 

Proof. Let F{x) = a;ln(x), x > 0. Then, F"{x) = l/x, i.e., F"{x) > for all x > 0. By 
Corollary [731 we get 



— / /(t)ln(/(f))At> -^ / /(t)Atlnf-^ / f{t)At) , 
-aJa b-aJa \b-aJa J 



and the result follows. D 

5.2 Applications to the Calculus of Variations 

We now show how the results obtained in Section lSTTJ can be applied to determine the minimum 
or maximum of certain problems of calculus of variations and optimal control on time scales. 

Theorem 80 (cf. [22] )• Let T be a time scale, a, b G T with a < b, and 99 : M — )• M &e a 
positive and continuous function. Consider the functional 



F[y{-)] 



ip{yit) + ht,it)y'^{t))dh\y^{t) 







At, a€M\{0, 1}, 



defined on all C^^ functions y : [a, 6]-^ — t- M satisfying y {t) > on [a,^]^; y(o) = 0, and 
y{b) = B. Define a function G{x) = J^ ip{s)ds, x > 0, and let G~^ denote its inverse. Let 

L (p(s)ds 
C= ^°/^ ^ . 5.7 

b — a 

(i) If a < or a > 1, then the minimum of F occurs when 

y{t) = G'\C{t-a)), te[a,%, 
and Fyain = {b — a)C°'. 
(ii) IfO<a<l, then the maximum of F occurs when 

y{t) = G-\C{t-a)), tG[a,6]T, 
and Fmax = {b - a)C°' . 
Remark 81. Since (^ is continuous and positive, G and G^^ are well defined. 
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Remark 82. In cases a = or a = 1 there is nothing to minimize or maximize, i.e., the 
problem of extremizing -F[y(-)] is trivial. Indeed, if a = 0, then F[y(-)] = b — a; if a = 1, then 
it follows from Theorem 1301 that 

{Goy)^{t)M 

= G{B) . 
In both cases -F is a constant and does not depend on function y. 
Proof of Theorem\8(A Suppose that a < or a > 1. Using Corollary 1751 we can write 



F[y{.)]>{b-a) 



l-a 



b ( rl 



ipiy{t) + hf,{t)y^{t))dh\y^{t)At 

= (6-a)i""(G(2/(6))-G(2/(a))r 



where the equality holds if and only if 

I / f{y{t) + hfi{t)y^{t))dh i y^{t) = c, for some c € M, t G [a, b]j. 

Using Theorem 1301 we arrive at 

iGoy)^{t)=C. 

A-integrating from a to t yields (note that y{a) =0 and G{0) = 0) 

G{y{t)) = c{t-a), 

from which we get 

y{t) = G-\c{t-a)). 

The value of c is obtained using the boundary condition y{b) = B: 

_ G{B) _ j^ ^{s)ds 



b — a b — a 

with C as in ()5.7p . Finally, in this case 

rb 



c, 



Fmin = / C"At ={b- a)C''. 
J a 



The proof of the second part of the theorem is done analogously using the second part of 
Corollary [751 □ 
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Remark 83. We note that the optimal solution found in the proof of the previous theorem 
satisfies y > 0. Indeed, 

y^{t) = {G^\C{t-a))f 

= [ iG~^y[C{t - a) + hfi{t)C]dh C 
Jo 

>0, 

because C > and (G~^)'(x) = (c-iix)) > ^ ^°^ ^^^ x >0. 

Theorem 84 (cf. [22]). Let if : [a,b]j -^M be a positive and rd-continuous function. Then, 
among all C^^ functions y : [a, b]j -^ M with y{a) = and y{b) = B, the functional 

F[y{-)]= f\{t)ey^('^At 

J a 

has minimum value Fmin = {b — a)e^ attained when 

y{t) = - I ln(v3(s))As + C{t-a), t£ [a, b]j, 

J a 



where 

b — a 
Proof. By Corollary [771 

F[y{-)]= t e^<^m+y^{t)^t 

J a 

with F{y{-)) = {b- a)e^[^" ^°^'^^*^^'^*+^] if and only if 

ln((^(t)) + y^{t) = c, for some c G M, t G [a, 6]| . (5.9) 

Integrating ()5.9p from a to t (note that y{a) = 0) gives 

y{t) = - I ln((^(s))As + c{t - a), t G [a, b]f . 

J a 

Using the boundary condition y{b) = B we have 

_ j'jn{y,{t))At + B _^ 

c — — o, 

b — a 

with C as in (j5.8p . A simple calculation shows that -Fmin = {b — a)e . D 

Remark 85. If we let T = M in the previous theorem we get |301 Theorem 3.4]. 
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Theorem 86 (cf. [22]). Let if : [a,b]j -^ IB. be a positive and rd-continuous function. Then, 
among all C^^-functions y : [a, b]j -^ M satisfying y > 0, y{a) = 0, and y{b) = B, with 

B+ f\(s)As 

{- ^^ ^ > y(t), tG[a,6]|, (5.10) 

— a 

the functional 

F[yi-)] = fw) + y^(t)] inb(t) + y^(t)]At 



has minimum value -Fmin = {b — a)Cln(C) attained when 

y{t) = C{t-a)- / ip{s)As, t G [a, bjj , 

J a 



where 



B+rip(s)As , , 

C = htZ\l — . 5.11 

b — a 



^[y(-)]> / k(t) + 2/"(t)]Atln(^ r[(^(t) + y^(t)]At 



Proof. By Corollary [79l 

u)]> fw)+y^m 

J a 

= ( ip{t)At + B\ In 

with F[y{-)] = (X,V(i)At + S) In (^^^^^^^^ if and only if 

ip{t) + y^{t) = c, for some c G M, t G [a, &]|. 
Upon integration from a to t (note that y(a) = 0), 

y{t) = c{t- a) - / V'(s)As, t G [a, 6]t- 



7>(t)At + S^ 



b — a 



Using the boundary condition y{b) = B, we have 



b — a 



s 



where C is as in (jS.lip . Note that with this choice of y we have, using (|5.10p . that y (t) = 
C - ip{t) > 0, t G [a, 6]|. It follows that Fmin = {b - a)Cln(C). D 

In order to close this subject we would like to point out that Theorem 3.6 in [30] is not 
true. This is due to the fact that the bound on the functional / considered in the proof is not 
constant. Let us quote the "theorem": 

47 



CHAPTER 5. INEQUALITIES APPLIED TO SOME CALCULUS OF VARIATIONS 

AND OPTIMAL CONTROL PROBLEMS 



Let 93 : M — 7- M be a positive and continuous function and a > 0. Then, among 

all C^ functions y : [0,a\ — > M satisfying y' > 0, y(0) = 0, and y{a) = A, the 

functional 

I = ln{{p{x)y' {x))dx 
Jo 

attains its maximum when 
where 

1 r 1 

C=- ^ds, (5.12) 

A Jo fis) 

and 

/max = -aln(C). (5.13) 

Now we give a counterexample to |30t Theorem 3.6]. Let a = A = 1, ip{x) = x + 1, and 
y{x) = X for all x G [0, 1]. Then, the hypotheses of the "theorem" are satisfied. Moreover, 

I[y{x)] = / \n{ip{x)y' {x))dx = [{x + l)(ln(2; + 1) - l)];;;io = 21n(2) - 1 « 0.386. 
JO 

According with (|5.12p and (|5.13|) the maximum of the functional / is given by /max = — In(C')) 
where 

C= / -—-ds. 
Jo f{s) 

A simple calculation shows that C = ln(2), hence /max = — ln(ln(2)) ~ 0.367. Therefore, 
I[y{x)] > /max, which proves our claim. 

5.2.1 An example 

Let us consider T = Z, a = 0, 6 = 5, -B = 25 and ip{t) = 2t + 1 in Theorem 1861 

Corollary 87 (cf. [22]). The functional 

4 
F[y{-)] = ^[{2t + 1) + {y{t + 1) - y(t))] ln[(2i + 1) + {y{t + 1) - y(t))], 
t=o 

defined for all y : [0, 5] n Z ^- M such that y{t + 1) > y{t) for all t S [0,4] fl T, attains its 

minimum when 

y{t) = lot -t^, te [0, 5] n Z, 

andFmin = 501n(10). 

Proof. First we note that max{c^(i) : t G [0, 4] n Z} = 9, hence 

b — a 5 

Observing that, when T = Z, (t^) = 2t + 1, we just have to invoke Theorem [86] to get the 
desired result. D 
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5.3 State of the Art 

The results of this chapter are aheady pubhshed [22]. Direct methods are an important 
subject to the calculus of variations theory and further research is in progress, extending 
Leitmann's direct method to time scales |73j . 
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Inequalities on Time Scales 



This chapter is devoted to the development of some integral inequalities as well to some of its 
applications. These inequalities will be used to prove existence of solution(s) to some dynamic 
equations and to estimate them, and this is shown in Sections 16.11 and 16.21 In Section [6.31 we 
prove Holder, Cauchy-Schwarz, Minkowski and Jensen's type inequalities in the more general 
setting of Oa-integrals. Finally, in Section [U31 we obtain some Gronwall-Bellman-Bihari type 
inequalities for functions depending on two time scales variables. 

Throughout we use the notations M.q = [0, oo) and M"*" = (0, oo). 



6.1 Gronwall's type inequalities 

We start by proving a lemma which is essential in the proofs of the next theorems. 

Lemma 88 (cf. [IS]). Leta,b G T, and consider a function r € C^^{[a,b]j,'U.'^) withr (t) > 
on [a,b]j. Suppose that a function g € C(M(J',MjJ') is positive and nondecreasing on M"*" and 
define, 

G{x)= r^ 

where x > 0, xq > if f^ -A < oo and x > 0, xq > if J^ -A = cx). Then, for each 
t € [a, b]j, we have 

G{r{t)) < G{r{a)) + £ ^(^Ar. (6.1) 
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Proof. Since g is positive and nondecreasing on (0,00), we have, successively, that 



r{t) < r{t) + /i^f(t)r^(t) , 
g{r{t)) < g{r{t) + hi^{t)r^{t)) , 



1 1 

< 



g{r{t) + hii{t)r^{t)) " g{r{t)) 
1 .. /-i 1 



(6.2) 



dh < / -—-—dh 



g{r{t)+hfi{t)r^{t)) ' Jq g{r{t)) g{r{t)) 

-dh\r-it)< "^"^'^ 



g{r{t) + h^,{t)r^it)) j ^ ^ " g(r(t)) ' 

for all t E [a, b]j and h G [0, 1]. By A-integrating the last inequality in (j6.'2|) from a to t and 
having in mind that Theorem [30] guarantees that 

(G o r)^(t) = I / G'{r{t) + h^i{t)r^{t))dh\ r^{t) 
-1 ^ ^ 

9(r(t) + Mt>^(i))^7'^^'^' 

we obtain the desired result [note that the case i = 6 if p{b) < 5 is also proved because of 
(US])]. D 



Theorem 89 (cf. [IS]). Let u{t) and f{t) be nonnegative rd-continuous functions in the time 
scales interval T* := [a, fe]T (ind TJ, respectively. Let k{t,s) be defined as in Theorem \29\ in 
such a way that k{t, s) and k^^{t, s) are nonnegative for every i, s G T* with s < t for which 
they are defined (it is assumed that k is not identically zero on T!^ x T^ ). Let $ G C(]R(J', RJ) 
be a nondecreasing, subadditive and submultiplicative function, such that <I>(u) > for u > 
and let W G C(M^,Mq') be a nondecreasing function such that for u > we have W{u) > 0. 
Assume that a{t) is a positive rd-continuous function and nondecreasing for t G T*. // 



uit)<a{t)+ I f{s)u{s)As+ I f{s)w(f A;(s,r)$(u(T))AT j As, 

Ja Ja \J a J 

for a<T<s<t<b, r, s, t G T*, then for all t G T* satisfying 

/•pit) ps 

*(C) + / Hpit), s)^ip{s))^{ / /(T)Ar)As G Dom(^-^) 

Ja J a 



(6.3) 



we have 



u{t) < p{t)a{t) 



+ p{t) j f{s)W ^~^U{0 + fk{s,rMpir))^(J^ fie)Ae]AT] 



As, (6.4) 
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where 



p(i) = l+ / fis)efit,ais))As, 

J a 



c = 



Pib) 



k{p{b),s)^{p{s)a{s))As, 
1 



and, as usual, ^"^ denotes the inverse of '^ . 



ds, X > 0, xo > . 



(6.5) 



(6.6) 



Remark 90. We are interested to study the situation when k is not identically zero on TJJ x T!J^ . 
That comprise the new cases, not considered previously in the literature. The case k{t, s) = 
was studied in jS] Theorem 3.1] and is not discussed here. 



Proof. Define the function z{t) in T* by 

rt 



z{t) = a{t)+ f{s)W[ A;(s,r)^(7/(T))ATJ As. 
Then, ()6.3p can be rewritten as 

u{t) < z{t) + / /(s)n(s)As. 

J a 

Clearly, z{t) is rd-continuous in i € T*. Using Theorem 1471 we get 

u{t) <z{t)+ I f{s)z{s)ef{t,a{s))As. 

J a 

Moreover, it is easy to see that z{t) is nondecreasing in t G T*. We get 

u{t) < z{t)p{t), 
where p{t) is defined by ()6.5p . Define 

v{t)= f k(t,s)^u{s))As, teT^ 

J a 

From (j6.8p . and taking into account the properties of ^, we get 



(6.7) 



(6.8) 



v{t) < / k{t,s)^ 



p{s) a{s) + / f{T)W{viT))AT 



As 



< 



< 



I k{t,s)^{p{s)a{s))As+ I k{t,s)<^(p{s) I /(T)W(w(r))Ar ) As 

Ja Ja \ Ja / 

fPib) 

/ k{p{b),s)<^{p{s)a{s))As 

+ / A:(t,s)$(p(s) / f{T)ATJ<l>{W{v{s))As 
C+ f kit,s)<^(p{s) I f{T)AT\^{W{v{s))As. 
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Define function r{t) on T^ by 

r{t) = C+ f k{t,s)<^(p{s) I /(r)Arj $(Vl/(t;(s))As. 

Since p and a are positive functions, we liave tliat <l>(a(s)p(s)) > for all s € T*. Since 
k^ > 0, we must have C, > 0, hence r(i) is a positive function on TJ. In addition, r(t) is 

2 

A-differentiable on TJ with 

+ / k^'{t,s)<^(p{s) I /(r)AT J $(VF(w(s))As 



k{a{t),t)^(^p{t) j^ fir) At 



(6.9) 



< $(W(r(t))) 

+ J fe^i (i, s)<^ (p{s) y" /(r) At") As 

Dividing both sides of inequality ()6.9p by $(VF(r(t))), we obtain 



^{W{r{t))) 



< 



j k{t,s)^(p{s) j' f{T)AT\As 



Let us consider the function ^ defined by (j6.6p . A-integrating this last inequality from a to 
t and using Lemma [551 we obtain 

^(r(t)) < -^{ria)) + / A;(t,s)$ U(s) / /(r)Arj As, 



5.10) 
D 



from which it follows that 

r{t) < ^-1 UiC) + J k{t, s)<^ip{s)Mj' f{T)Ar)As^ , t G T^ (l 

Combining ()6.10p . ()6.8p and ()6.7p . we obtain the desired inequality (|6.4p . 

If we let T = M in Theorem [89} we get |35l Theorem 2.1]. If in turn we consider T 
then we obtain the following result: 



Corollary 91 (cf. |46] ) . Let u{t) and f{t) he nonnegative functions in the time scales interval 
T* := [a, 6]z and [a,b — l]z, respectively. Letk(t,s) be defined as in Theorem\2^in such a way 
that k{t, s) and k '^{t, s) = k{a{t), s) — k{t, s) are nonnegative for every i, s € T* with s < t for 
which they are defined (it is assumed that k is not identically zero on [a, b — 1]t» x [a, b — 2]j_^). 
Let <1> € C(M(]',M^) be a nondecreasing, subadditive and submultiplicative function such that 
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^{u) > for n > and let W S C(M(|',M(j') be a nondecreasing function such that for u > 
we have W{u) > 0. Assume that a{t) is a positive and nondecreasing function for t G T*. // 

t-i t-i /s-l \ 

u{t) < a{t) + J2 fi^H^) + E /(*)^ E ^(^' r)Hn{T)) , 

s=a s=a \T=a / 

for a<T<s<t<b, T,s,t£T^, then for all t G T* satisfying 

t-2 s-l 

^(C) + Y.^{t- 1, s)^{p{s))^Y. /(^)) ^ Dom(^-i) 



we have 



t-i 



t{t)<p{t)\a{t) + Y,f{s)W 



T-l 



^^' h'(^) + E ^(^' ^)^(Kr))$(E / w) 



6»=a 



where 



t-i 



p{t) = l + Y,f{s)ef{t,s + l) 



6-1 



^ = ^fc(6-l,s)$(p(s)a(s)), 



ond ^^""^ is i/ie inverse of '^ . 



s=a 
1 



,■0 ^{W{s)) 



ds, X > 0, rro > , 



For the particular case T = M, Theorem [92] generahzes the result obtamed by Oguntuase 
in Uni Theorems 2.3 and 2.9]. 

Theorem 92 (cf. [46j ) . Suppose that u{t) is a nonnegative rd-continuous function in the 
time scales interval T* = [a, 6]^ and that h{t), f{t) are nonnegative rd-continuous functions 
in the time scales interval T^. Assume that b{t) is a nonnegative rd-continuous function and 
not identically zero on TJ . Let ^{u), W{u), and a{t) be as defined in Theorem \89[. If 

u{t) < a{t) + I f{s)u{s)As + I f{s)h{s)W ( I 6(T)$(n(r))Ar ) As, 

Ja Ja \J a / 

for a < T < s < t <b, T,s,t (^ T^,, then for all t G T^, satisfying 

-Pit) 



rP(t) t-T 

*(6 + / fc(T)$(p(r))$( / f{e)h{e)A9)AT G Dom(*'i) 

Ja J a 



we have 



u{t) < pit)a{t) 



+ pit) / f{s)h{s)W 



^~' ^(0 + / b{r)Hp{r))H / f{e)h{e)A9)AT 



As. 
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where p{t) is defined by \6.5\) . ^ is defined by Ii6.6\) . and 

rpib) 
^= b{s)^{p{s)a{s))As. 

J a 

Proof, similar to the proof of Theorem [89l D 

For the remaining of this section, we use the fohowing class of S fmictions. 

Definition 93 {S function). A nondecreasing continuous function g : W^ — > M^ is said to 
belong to class S if it satisfies the following conditions: 

1. g{x) is positive for x > 0; 

2. {1/ z)g{x) < g{x/z) for x > and z > 1. 

Remark 94. For a brief discussion about this class of S functions, the reader is invited to 
consult [191 Section 4]. 

Theorem 95 (cf. [l6]). Let u{t), f{t), k{t,s), <& and W be as defined in Theorem\8M and 
assume that g & S. Suppose that a{t) is a positive, rd-continuous and nondecreasing function. 
If 

u{t) < a{t) + f f{s)g{u{s))As + f f{s)W ( f k{s, t)'^{u{t))At) As, (6.11) 

Ja Ja \J a J 

for a<T<s<t<b, r, s, t € T^,, then for all t € T^, satisfying 

Gil) + / /(r)AT G Dom(G-i) 

J a 

r-piyt) r-T 

*(C) + / k{p{t),T)^{q{T))^{ / f{e)M)AT G Dom(M/-l), 

J a J a 



and 

rpit) 



we have 



u{t) < q{t) max{a(t), 1} 

+ q{t) [ f{s)wh-'U{C)+ rk{s,TMq{T)Mr f{0)Ae)A, 

Ja L V Ja Ja 



As. 



where ^ is defined by i6.6\} . 



f^ ds 
G{x) = / — -, X > 0, 5 > 0, 
Js 9{s) 

r-t 



q{t) = G-' (^G(l) + ^ /(t)At) , (6.12) 

fP{b) 

C= / k{p{b),s)<5{q{s)max{a{s),l})As, 

J a 

and G^^ is the inverse function of G. 
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Proof. Define the function 

z(t) =max{a(t),l}+ / f{s)w(f A:(s, t)$(u(t))At J As. 
Then, from ()6.1ip we have that 

uit) < z{t) + f fis)giuis))As. 

•la 

Clearly, z{t) > 1 is rd-continuous and nondecreasing. Since g & S, we have 



z{t) 



or 



xit)<l+ f f{s)g{x{s))As, (6.13) 

J a 

with x{t) = u{t)/z{t). If we define v{t) as the right-hand side of inequality ()6.13p . we have 
that v{a) = 1, 

v^{t) = f{t)g{x{t)), 

and since g is nondecreasing, 

v^{t) < f{t)g{v{t)), 

or 



Being the case that v (t) > 0, A-integrating (j6.14p from a to t and applying Lemma [881 we 
obtain 

Givit)) < G(l) + / /(r)AT, 

J a 

which implies that 

v{t)<G''(^G{l) + J^ /WAt^ 

We have just proved that x{t) < q{t), which is equivalent to 

u{t) < q{t)z{t). 

Following the same arguments as in the proof of Theorem 1891 we obtain the desired inequality. 

D 

If we consider the time scale T = /iZ = {hk : k G Z}, where h > 0, then we obtain the 
following result. 
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Corollary 96 (cf. [36]). Let a,b £ hZ, h > 0. Suppose that u{t), f{t), k{t,s), $ and W 
are as defined in Theorem [221 and assume that g £ S. Suppose that a{t) is a positive and 
nondecreasing function. If 

u{t)<a{t)+ Yl f{s)g{u{s))h+ J] f{s)wi ^ k{s,T)^u{T))h\ h, 

«S[a,t)T, sG[a,t)T, \T"e[a,s)T* / 

for a < T < s < t <b, T,s,t gT^: = [a, b]hz, then for all t G T^, satisfying 



and 



we have 



G(l)+ Yl f(.r)heBom{G~') 

re[a,t)T« 



*(C)+ Yl k{t-h,T)^q{T))^i Y f{0)h]hGDom{^-^) 

Te[a,t-/i)T, \06[a,T)T, 



u{t) < q{t){ max{a(i),l} 



+ Y /(«)^ 



TG[a,s)T^ \6»e[a,T)Tr, 



where ^ is defined by i6.6\} 



G{x) = / — -, X > 0, 6>0, 
Js 9[s) 



q{t) = G-nG{l) + Y /(^)^ 
C= Y^ k{b- h,s)^{q{s)max{a{s),l})h, 

se[a,b-h)T:t 

and G~^ is the inverse function of G. 

Theorem 97 (cf. [56]). Letu{t), f{t), b{t), h{t), $ andW be as defined in Theorem [M and 
assume that g £ S. Suppose that a{t) is a positive, rd-continuous and nondecreasing function. 
If 

n(t)<a(t)+ f f{s)g{uis))As+ f f{s)h{s)W ( T b{TMu{T))AA As, 

Ja Ja \J a J 

for a < T < s < t <b, T,s,t £ T^,, then for all t € T^, satisfying 

*(0+/ Kr)^(.Q{r))^(. f{9)h{9)Ae)AT£Bom{^~^), 
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we have 

u{t) < q(t) inax{a(t), 1} 

+ q{t) f f{s)h{s)W 



\ J a J a 



As, 



where ^ is defined by 116. 6\) . q{t) is defined by Ii6.1^) and 

fPib) 

^= 6(s)$(g(s)max{a(s),l})As, 

J a 

Proof, similar to the proof of Theorem [95l D 

We now make use of Theorem [22] to the quahtative analysis of a nonlinear dynamic 
equation. Let a, 6 € T and consider the initial value problem 

u^{t) = F (t,u{t), K{t,u{s))As\ , t£T'^, u{a)=Ua, (6.15) 

where T^ = [a,b]T, u £ Cl^{T^), F G Crd(T* x M x M,M) and K G Crd(T* x R,R). 

In what follows, we shall assume that the IVP (|6.15p has a unique solution, which we 
denote by u^{t). 

Theorem 98 (cf. HB]). Assume that the functions F and K in \6.15\) satisfy the conditions 

\K{t,u)\<h{t)^{\u\), (6.16) 

\F{t,u,v)\ < \u\ + \v\, (6.17) 

where h and $ are as defined in Theorem \92[ Then, for t G T^, such that 

/■P(t) fr 

^(0+/ ^ip{r))^{ /i(e)A0)AT G Dom(^-^), 

J a J a 

we have the estimate 

\u,{t)\ < p{t)!.\ua\ + I /i(s)^-i U{i) + r <^{p{t))<^{ r h{9)Ae)AT] As\ , (6.18) 



where 



rt 
p{t) = l+ I ei{t,(7{s))As, 



rpy.0) 
i= ^{p{s)\ua\)As, 

J a 

^{x) = I -^rrrds, x > 0, xq > 0. 
JxQ ^(s) 
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Proof. Let n*(i) be the solution of the IVP (|6.15p . Then, we have 

u4t)=Ua+ f(s,u^{s), K{s,u^{t))At)] As. (6.19) 

Using ([636]) and (l6T71l in (fOOD . we have 

^a\+ (\Ms)\+ \K{s,u^{t))\Atj As 

I (\u,{s)\ + h{s) I $(|n,(r))|)Arj As. (6.20) 



ft 

\uA < \Ua 



(■t 
< \Ua\ + 



A suitable application of Theorem [U2] to (|6.2U|) . with a{t) = \ua\, f{t) = b{t) = 1 and 
W{u) = u, yields (fHJHD . D 

6.2 Some more integral inequalities and applications 

In this section we shall be concern with the existence of solutions of the following integrody- 
naniic equation 

x^{t)=F(t,x{t), K[t,T,x{T)]AT], x{a) = A, tG[a,6]f, (6.21) 

where a,b gT, K : [a, b]j x [a, b]j xR ^ M. and F : [a, 6]| x M^ — )• M are continuous functions. 

As is well known integrodifferential equations and their discrete analogues find many 
applications in various mathematical problems. Moreover, it appears to be advantageous to 
model certain processes by employing a suitable combination of both differential equations 
and difference equations at different stages in the process under consideration (see [67j and 
references therein for more details). 

Some results concerning existence of a solution to some particular cases of the integrody- 
namic equation in ()6.2ip were obtained in |63] [95] . Here we will make use of the well-known 
in the literature Topological Transversality Theorem to prove the existence of a solution to 
the above mentioned equation. 

We first introduce some basic definitions and results on fixed point theory. Let ;S be a 
Banach space and C C B he convex. By a pair {X, A) in C is meant an arbitrary subset X 
of C and an A C X closed in X. We call a homotopy H : X x [0, 1] -^ Y compact if it is 
a compact map. If X C Y, the homotopy H is called fixed point free on A C X if for each 
A G [0, 1], the map H\A x {A} : A ^ Y has no fixed point. We denote by Ca{X, C) the set of 
all compact maps F : X —^ C such that the restriction F\A : A ^ C is fixed point free. 

Two maps F,G G CAiX,C) are called homotopic, written F ~ G in Ca{X,C), provided 
there is a compact homotopy Hx : X ^>- C (AG [0, 1]) that is fixed point free on A and such 
that Ho = F and Hi = G. 

60 



6.2. SOME MORE INTEGRAL INEQUALITIES AND APPLICATIONS 

Definition 99. Let {X,A) be a pair in a convex C C B. A map F G CAiX,C) is called 
essential provided every G S CAi^, C) such that F\A = G\A has a fixed point. 

Theorem 100 (Topological Transversality [SI]). Let {X,A) be a pair in a convex C C B, 
and let F, G be maps in Cyi(X, C) such that F c:^ G in Ca{X, C). Then, F is essential if and 
only if G is essential. 

The next theorem is very useful to the application of the Topological Transversality The- 
orem. Its proof can be found in 



Theorem 101. Let U be an open subset of a convex set C C B, and let {U,dU) be the pair 
consisting of the closure of U in G and the boundary of U in G. Then, for any uq G U , the 
constant map F\U = uq is essential in Cqu{U, G). 



If more details are needed about fixed point theory we recommend the books 
Here, as usual, G{[a., 6]t, M) [sometimes we will only write G[a., 6]t] denotes the set of all real 
valued continuous functions on [a, b\j and G^{[a, 6]^, M) the set of all A-differentiable functions 
whose derivative is continuous on [a,6]|, and we equip the spaces C([a, 6]Tr,K), G^{[a,b]f,'M) 
with the norms ||n||o = supjg[^_b]^ \u{t)\, \\u\\^ = suptg[„^b]^ \u{t)\ + supjg[„^b]«; \u^{t)\, respec- 
tively. 

To make use of the Topological Transversality Theorem it is essential to obtain a priori 
bounds on the possible solutions of the equation in study. To achieve that, the next lemma 
is indispensable. 

Lemma 102. Let f,g € C([a, 5]^, Rq") and u,p G C([a, 6]Tr,IK|^), being p{t) positive and 
nondecreasing on [a,b]j, k G G{[a,b]j x [a, 5]| ,R(J^), and c G M,|. Moreover, let w,w G 
C(Mq',M^) be nondecreasing functions with {w,w}{x) > for x > 0. Let B,D > and 
M = [b — a)D. Define a function G by. 



G{x) = / — — , ... ds, X > 0, 

Ji w(Taa.yi{s,Mw[s)\) 

and assume that lim^^^^oo G(x) = oo. If, for all t G [a,6]Tr, the following inequality holds 

u{t)<p{t) + B w( max <u{s) + c,D w{u{t) + c)At>] As, 

then, 

u{t) < G'^ [G{p{t) + c) + B{t -a)]-c, t G [a, b]j. 

Proof. We start by noting that the result trivially holds for t = a. Consider an arbitrary 
number to ^ (a, ^]t and define a positive and nondecreasing function z{t) on [a, toll by 

z{t) = p{to) + B w ( max < u{s) + c, D w{u{t) + c) At I j As. 
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Then, z{a) = p(to)) for all t S [a, toll we have u{t) < z{t), and 

z^{t) = Bw i max \ u{t) + c, D w{u{t) + c)Ar 

< Bw (max {u{t) + c, D{b - a)w{z{t) + c)}) 
= Bw (max {u{t) + c, Mw{z{t) + c)}) , 

for all t S [a, tolif- Hence, 

z^{t) 



<B, 



w (max {u{t) + c, Mw{z{t) + c)}) 
and, after integrating from a to t and with the help of Lemma [551 we get 

G{z{t) + c) < G{z{a) + c) + B{t-a), te [a, to]T- 

In view of the hypothesis on function G, we may write 

z{t) < G~^ [G{z{a) + c) + B{t - a)] - c. 

Therefore, 

uit) < G'^ [G{pito) + c) + B{t - a)] - c, 

for all t G [a, toll- Setting t = to in the above inequality and having in mind that to is 
arbitrary, we conclude the proof. D 

Remark 103. We note that by item 1 of Theorem [7] there is an inclusion of G^[a,b]j into 
G[a,b]j. 

Theorem 104. The embedding j : C^[a, &]t — )• C[a, &]t is completely continuous. 

Proof. Let i? be a bounded set in C^[a, &]t. Then, there exists M > such that ||x||i < M 
for all X (z B. By the definition of the norm || • |ji, we have that sup^gr^ u 1^(^)1 ^ -^) hence 
ll^^llo ^ M, i.e., B is bounded in C[a, b]j. Let now e > be given and put 6 = -^. Then, it is 
easily seen that, for x ^ B we have |x (t)| < M for all t € [a, b]j. For arbitrary ti, t2 £ [o-, bjj 
with ti 7^ ^2) we use Theorem [31] to get 

\xiti) - X{t2)\ ^^ 



\tl-t2\ 

i.e., for \ti — t2\ < S we have \x{ti) — x{t2)\ < e and this proves that B is equicontinuous on 
C[a,b]j (the case ti = t2 is obvious). By the Arzela-Ascoli Theorem, B is relatively compact 
and therefore j is completely continuous. D 
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Theorem 105. Assume that the functions F and K introduced in ()6.2ip satisfy 

\F{t,x,y)\ <S«;(max{|x|,|y|}) + C, (6.22) 

\K{t,s,x)\<Dw{\x\), (6.23) 

where B, C and D are positive constants and w,w are defined as in Lemma \102[ 
Then, the integrodynamic equation in ()6.2ip has a solution x € C^([a, 6]^,!^). 

Proof. We start considering the following auxiliary problem, 

y^{t) = FUy(t)+A, j K[t,s,y{s)+A]As\, y(a) = 0, tG[a,6]|, (6.24) 

and showing that it has a solution. To do this, we first establish a priori bounds to the 
(possible) solutions of the family of problems 

y^{t) = \FUy{t) + A, j K%s,y{s) + A]^^, y(a) = 0, tG[a,6]|, (6.25) 

independently of A € [0, 1] . 

Integrating both sides of the equation in ()6.25p on [ajt]-!-, we obtain, 

y{t) = xJfU y{s) + A,j K[s, r, 2/(r) + A]AtJ As, 

for all t G [a, b]f. Applying the modulus function to both sides of the last equality we obtain, 
using the triangle inequality and inequalities ()6.22p and ()6.23p . 

\yit)\ <p{t) + B I u;fmax||y(s) + A|,D / w{\y{T) + c\)AtX\ As, 

with p{t) = C{t — a). By Lemma [1021 (with u{t) = |y(t)|, c = \A\), we deduce that (note also 
that \y{.) + A\<\y{-)\ + \A\) 

\y{t)\ < G-^ [G{p{t) + 1^1) + B{t - a)] - \A\, (6.26) 

for all t G [a, 6]t. Denote the right-hand side of ()6.26p by R{t). 

Let B = C^[a,h]'f be the Banach space equipped with the norm || • ||i and define a set 
C = {n G C^la, h\f : u{a) = 0}. Observe that C is a convex subset of B. 

Define the linear operator L : C ^ C[a, 6]t by Lu = u . It is clearly bijective with a 
continuous inverse L~^ : C[a, b]j — t- C. Let Mq and Mi be defined by 



Mo = Rib), 
Ml = sup 

t6[a,fe]T,|y|<Mo 



F[t,y, I K[t,s,y]As 
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Consider the family of maps Tx ■ C[a, b]j -^ C[a, b]j, < A < 1, defined by 



iJ'xym 



XF U y{t) +aJ K[t, s, y{s) + A]As\ , 



and the completely continuous embedding j : C ^ C[a, b]r (it is easily seen that the restriction 
here used of j of Theorem I1U4I is also completely continuous) . Let us consider the set Y = 
{2/ S C : ||y||i < r} with r = 1 + Mq + Mi. Then we can define a homotopy Hx : Y ^ C 
by Hx = L^^j^xj- Since L~^ and Fx are continuous and j is completely continuous, H \s a. 
compact homotopy. Moreover, it is fixed point free on the boundary of Y . Since Hq is the 
zero map, it is essential. Because Hq ~ Hi, Theorem 1 100 1 implies that Hi is also essential. 
In particular. Hi has a fixed point which is a solution of ()6.24p . To finish the proof, let 
y G C^la, h\j be a solution of (|6.24p and define the function x{t) = y{t) + A, t ^ [o, b]j. Then, 
it is easily seen that x(a) = A and 

x^it)=F(t,x{t),j K[t,s,xis)]AsY tG[a,6]|, 

i.e., X G C"^[a, 6]ir is a solution of ()6.21|) . D 

Theorem 11051 was inspired by the work of Adrian Constantin in [S2] . It provides a non- 
trivial generalization to the one presented there, and as a consequence. Theorem 11051 seems 
to be new even if we let the time scale to be T = M. Indeed, let 

/ (x + l)ln(x + l) ifx>0; 

F{t,x,y) = \y\ + i -r n 

[0 if x < 0. 

We win show that, for ah /, /i G C([a, 6]|,IR^), \F{t,x,y)\ > f{t)\x\ + \y\ + h{t) for some 
{t,x,y) G [a,b]j x R^ [in [35] the assumption on F was that \F{t,x,y)\ < /(t)|x| + \y\ + h{t)]. 
For that, suppose the contrary, i.e., assume that there exist f,h ^ C([a, 5]^, Mq") such that 
\F{t,x,y)\ < f{t)\x\ + \y\ + h{t) for all {t,x,y) G [a,&]| x M^. In particular, for arbitrary 
X > 0, we have that, 

\y\ + {x + l) ln(x + 1) = \F{t, x,y)\ < f{t)x + \y\ + h{t), 

which is equivalent to 

(x + 1) ln(x + 1) < f{t)x + h{t), 

or 

x[ln(x + 1) - f{t)] + ln(x + 1) < h{t). (6.27) 

Now, we fix t G [a, b]j and let x ^' 00 in both sides of ()6.27p getting a contradiction. 

Note that the function w{x) = (x + 1) ln(x + 1), x > 0, is nondecreasing and is such that 
f?^ -rry-^ = oo . MorBOVcr, for F as above, we have 

Jl w(s) ' ' 

\F{t,x,y)\<\y\ + w{\x\) 

< max{|x|, \y\} + u'(max{|x|, \y\}). 
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If we define w{x) = max{x, w{x)}, x > 0, then w{x) < x + w{x) and (see [31]) /i°° ^;;;t^^^ = ^^ 
wliich shows that we can apply our result to such a function F. 

Corollary 106. Assume that p € C{[a,b]j,M). Then, the initial value problem 

x^{t)=p{t)x{t), x{a) = A, (6.28) 

has a unique solution x € C'^([a, fe]T,IR). 

Proof. Define F{t,x,y) = p{t)x, for (t, x) € [a,b]j x M. Then, \F{t,x,y)\ = \p{t)\\x\ and the 
existence part follows by Theorem 11051 

Suppose now that x,y € C^{[a,b]j,'M.) both satisfy (|6.28p . An integration on [a,t]j yields, 

x{t) = A+ p{s)x{s)As, 

J a 

y{t)=A+ I p{s)y{s)As. 
J a 

Hence, 

\x{t)-y{t)\< I \p{s)\\x{s)-y{s)\As. 



Using Gronwall's inequality (cf. Theorem B7|) [note that \p{s)\ S TV^] it follows that \x{t) — 
y{i)\ ^ and finally that x{t) = y[t) for all t G [a, b]f. Hence, the solution is unique. D 

Remark 107. S. Hilger in his seminal work [58] proved (not only but also) the existence of a 
unique solution to equation (|6.28p with p G Crd([a, &]|,K) being regressive. In this paper we 
are requiring p{t) to be continuous but not a regressive function. 

In view of the previous remark it is interesting to think of what happens when the function 
p is not regressive. This is shown in the following result. 

Proposition 108. Suppose that p € C([a, &]^,]R) is not regressive. Then, there exists io ^ 
[a, 6]| such that the solution of (j6.28p satisfies x{t) = for all t G [a{tQ),b]'f. 

Proof. Since p is not regressive, there exists a point to G V^-, ^It such that 1 + p{to)p{to) = 0. 
This immediately implies that /i(to) / and p(to) 7^ 0. Let x G C^([a, 5]t,K) be the unique 
solution of (|6.28p . Using formula (|1.2p it is easy to derive from (|6.28p that x^(t)[l+/i(t)p(t)] = 
p{t)x{a{t)) and in particular that x {to)[l + fi{to)p{to)] = p{to)x{a{to)). It follows that 
x(cr(to)) = 0. Now note that, since to is right-scattered, we must have o"(io) 7^ a. Moreover, 
x{t) = satisfies x'^{t) = p{t)x{t) for all t € [a{to),b]f. The uniqueness of the solution 
completes the proof. D 
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6.2.1 Some particular integrodynamic equations 

We now use Theorem I1U5I to prove existence of solution to an integral equation describing 
some physical phenomena. We emphasize that, for each time scale we get an equation, i.e., 
we can construct various models in order to (hopefully) better describe the phenomena in 
question. 

Theorem 109. Let 0,6 G T with 6 > and L € C^([0, 6]Tr,Ko'). Moreover, assume that 
a function M[t,s) G C([0,6]Tr x [0,6]|,M^) has continuous partial IS.- derivative with respect 
to its first variable [denoted by M *{t,s)] and M{a{t),t) is continuous for all t € [0,6]|. 
Suppose that the function w € C(M(|,]R(j') is nondecreasing, such that w{x) > for all x > 
and J-|°° :^r^ds = oo. Then, the integral equation, 

x{t) = L{t)+ M{t,s)w{x{s))As, tG[0,b]j, (6.29) 

Jo 

has a solution x G C"^([0, 6]T,IR(|)• 
Proo/. Let us define 

Fit, X, y) = L^{t) + M{a(t),t)w(\x\) + y, [t, x, y) G [0, b]^ x M^, 
K{t,s,x) = M^'{t,s)w{\x\), {t,s,x) G [0,6]| x [0,6]|' x M. 

Then, we have 

\F{t,x,y)\ < \L^{t)\+M{a{t),t)iv{\x\) + \y\, 
\Kit,s,x)\<\M^'it,s)\iv{\x\). 

Using Theorem 11051 with 



B > max M(a(t),t), 

te[o,fe]« 

C> max \L^{t)\, 

t&[0,b]j 

D> max \M^'{t,s)\, 

(t,s)e[0,6]«x[0,b]«^ 

'w{x) = id{x) and w{x) = max{x,u;(a;)} we conclude that the equation 

x^{t) = L^{t) + M{a{t),t)iD{\x{t)\) + I M^'{t,s)w{\x{s)\)As, 

Jo 

with initial value x{0) = L{0) has a solution x G C^([0, 6]t,IK). Now, an integration on [0,4]^ 

gives, using Lemma [221 

r-t 

xC] "['] 

'0 

By the assumptions on functions L, M and w we conclude that x{t) > for all t G [0,6]t, 
hence x is a nonnegative solution of ()6.29p . D 
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;(t) = L{t) + / M{t, s)w{\x{s)\)As, t G [0, 5]t. 
Jo 



6.2. SOME MORE INTEGRAL INEQUALITIES AND APPLICATIONS 

Corollary 110. If in equation ()6.29p . L{t) > on [0, 6]^ and w = w is continuously differ- 
entiable on (0,cx3), then the solution obtained by Theorem \109\ is unique. 

Proof. Assume that there exist two positive solutions x,y on [0,6]^ to ()6.29p . Then, 

x{t) - yit) = I M{t, s)[w{x{s)) - w{yis))]As, t G [0, b]j, 
Jo 

hence, 

\x{t) - y{t)\ < I M{t, s)\wix{s)) - w{y{s))\As, t G [0, b]j. 
Jo 

Let us now define the following numbers: 

71 = min {x{s),y{s)}, 

se[0,b]-r 

^2 = max {x{s),y{s)}, 

se[o,b]t 

V = max w' [x), 

a;e[7i.72] 

/i = max M[t, s). 

(t,s)G[0,fe]TX[0,6]^ 

The mean value theorem guarantees that 

\w{x{t)) - w{y{t))\ < u\{x{t) - y{t))\, t G [0, 6]t, 
hence, 



\x{t)-y{t)\< ^iu\x{s) - y{s)\As, tG[0,6]T. 
Jo 

Using Theorems?! we conclude that \x{t)—y{t)\ < on [0, b]f and this implies that x{t) = y{t) 
on [0,&]t. □ 

Now we want to mention some particular cases of equation (|6.29|) . Define on Mq" the 
function w{x) = x"^ , < r < 1. Then, the equation 

x{t) = L{t) + I M{t,s)[x{s)YAs, t G [0,6]t, (6.30) 

Jo 

has a unique positive solution if L{t) > for all t G [0, 6]t- This type of equation appears 
in many applications such as nonlinear diffusion, cellular mass change dynamics, or studies 
concerning the shape of liquid drops [72j • Some results concerning existence and uniqueness 
of solutions to (|6.30p [the time scale being T = R] were previously obtained (consult [79] and 
references therein). 

If we define u[t) = [x[t)Y , it follows from ()6.30p that, 

r-t 

[u{t)]r = L{t) + M{t,.s)u{s)As, tG[0,6]T. (6.31) 

Jo 
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It is easily seen that u{t) is the unique solution of ()6.3ip and, if we let r = | and M{t, s) = 
K{t - s) for K G C^(M([,M([), it follows that 

[u{t)f = L{t) + K{t - s)u{s)As, te[0,b]T. (6.32) 

Jo 

This equation appears in the mathematical theory of the infiltration of a fiuid from a cylin- 
drical reservoir into an isotropic homogeneous porous medium |8H [82] . Some existence and 
uniqueness theorems regarding solutions of (j6.32p were obtained in [32l [33l [Ml [SD [82] . 

Remark 111. Corollary II 101 proves uniqueness of a solution to the integral equation ()6.32p 
under different assumptions on the function L than those in [321 [331 [311 [SU [52] (obviously 
considering T = M). For example, (to prove uniqueness) in [33], the author considered 
L G C^(M+,M+) such that L'(0) / 0. Therefore, if we, e.g., let L{t) = t^ + l, t G [0,6], we see 
that L'{0) = 0, hence we cannot use the results obtained in [33]. If we let L{t) = t, Corollary 
IllOl cannot be applied. 

We close this section with an example of what the integrodynamic equation in (|6.2ip could 
be on a particular time scale different from M, namely, in T = hZ = {hk : A: E Z} for some 
h > 0. Let a = A = and b = hm for some m G N. Then, on this time scale, (|6.2ip becomes 

^i^±l^lz^ = FUx{t),Y,hK[t,kh,x{kh)]\, x(0)=0, t€[a,b]l^. 
\ fc=o / 

6.3 Oa-integral inequalities 

In this section we propose to prove inequalities of Jensen's, Minkowski's and Holder's type 
using the ()a~ integral (cf. definition in (j6.34p below). 

Based on the A and V dynamic derivatives, a combined dynamic derivative, so called Oa- 
derivative, was introduced as a linear combination of A and V dynamic derivatives on time 
scales [90], i.e, for < a < 1, 

/^'^(i) = a/^(t) + (l-Q)/^(t), tGi;:, (6.33) 

provided / and / exist. The Oo-derivative reduces to the A-derivative for a = 1 and to 
the V-derivative for q = 0. On the other hand, it represents a "weighted dynamic derivative" 
on any uniformly discrete time scale when a = g- A motivation for study such derivatives is 
for designing more balanced adaptive algorithms on nonuniform grids with reduced spuriosity 

m- 

Let a, i G T, and /i : T — )■ M. Then, the Oa-integral of h from a to i is defined by 

I h{T)<)aT = a h{T)AT + {l-a) h{T)\/T, 0<a<l, (6.34) 

Ja Ja Ja 
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provided that there exist A and V integrals of h on T. It is clear that the diamond-a integral 
of h exists when /i is a continuous function. We refer the reader to [701 1871 190] for an account 
of the calculus associated with the Oa-derivatives and integrals. 

The next lemma provides a straightforward but useful result for what follows. 

Lemma 112 (cf. [H]). Assume that f and g are continuous functions on [a, 6]t- If f{t) ^ oit) 
for all t G [a,b]f, then !lf{t)()^t < f^g{t)Oat. 



Proof. Let f{t) and g{t) be continuous functions on [a, 6]t. Using (|1.5p and the analogous in- 
equality for the V-integral we conclude that f^ f{t)^t < f^ g{t)At and f^ /(t)Vt < /^ g(t)Vt. 
Since a G [0, 1] the result follows. D 



We now obtain Jensen's type Oa-integral inequalities. 

Theorem 113 (cf. [H]). Let T be a time scale, a, 6 € T with a < b, and c, d G M. // 

g G C{[a,b]f,{c,d)) and f G C((c, (i),M) is convex, then 

\ b — a J b — a ' 

Remark 114. In the particular case a = 1, inequality ()6.35p reduces to that of Theorem 1721 
If T = M, then Theorem 11131 gives the classical Jensen inequality, i.e., Theorem 1441 However, 
if T = Z and f{x) = — ln(x), then one gets the well-known arithmetic-mean geometric-mean 
inequality (|5.39|) . 

Proof. Since / is convex we have 

Using now Jensen's inequality on time scales (cf. Theorem I72p. we get 

y b-a J b-a Ja b-a J^ 

= 6^("/ f{gis))As + {l-a)J^ /(5(s))Vs) 
f{gis))0as. 

D 



1 rb 



b-a.,^ 



Now we give an extended Jensen's inequality (cf. [TT]). 
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Theorem 115 (Generalized Jensen's inequality). Let T be a time scale, a, 6 € T with a < b, 
c, d eR, g e C{[a, 6]t, (c, d)) and h G C{[a, b]j,R) with 

b 
\Ks)\OaS >0. 



/// G C((c, d),M) is convex, then 

. ( la\HsMs)Oas \ ^ j'jh{s)\f{g{s))O^S 
V Ia\Hs)\OaS J - l'jh{s)\OaS 

Remark 116. In the particular case h = 1, Theorem 11151 reduces to Theorem 11131 

Remark 117. Similar result to Theorem II 151 holds if one changes the condition "/ is convex" 
to "/ is concave", by replacing the inequality sign "<" in (j6.36p by ">". 

Proof. Since / is convex, it follows, for example from |5H Exercise 3.42C], that for t G (c, d) 
there exists a^ G M such that 

atix -t)< f{x) - f{t) for all x G (c, d). (6.37) 

Setting 

^_ ta\Ks)\9{s)0o.S 
ta\Ks)\O.S ' 

then using ()6.37p and Lemma 11121 we get 

Hs)\f{9{s))0aS - / \h{s)\OaS f , 

= f \h{s)\f{g{s))Oo^s-( f \his)\Oas)f{t)= f \h{s)\{f{g{s))-f{t))0as 
>at(f \h{s)\{g{s)-t)\0as = at( f |/i(s)|g(s)OaS - t /" \h{s)\OaS 

=at(f \h{s)\g{s)0as- f \h{s)\g{s)^^s^ = . 

This leads to the desired inequality. D 

Now we sate two (well-known in the literature) corollaries of the previous theorem. 

Corollary 118. (T = M) Let g,h : [a,b] — )• M 6e continuous functions with g{[a,b]) C (c, d) 
and J \h{x)\dx > 0. If f & C((c, d),M) is convex, then 

f f^\h{xMx)dx \ ^ f^\h{x)\f{g{x))dx 
[ J^\h{x)\dx J - J^\h{x)\dx 
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Corollary 119. (T = Z) Given a convex function f, we have for any xi, . . . ,Xn S M and 

ci, . . . ,c„ G M with Ylk=i \^k\ > 0/ 



/ [%^^) < ^^='J"i^^"'^ ■ (6.38) 



Particular Cases 



(i) Let g(t) > on [a, feji and /(t) = t^ on (0, +oo). One can see that / is convex on 
(0, +oo) for /3 < or /3 > 1, and / is concave on (0, +oo) for /3 € (0, 1). Therefore, 



(ii) Let g{t) > on [a, 6]^ and fit) = ln(t) on (0,+oo). One can also see that / is concave 
on (0, +00). It follows that 



^^ fa\h{sMs)OaA ^ fjh{s)\ln{g{s))^^s 



(iii) Let h = 1, then 



^^f i:9{s)Ws \f^Hg{s))0, 



b — a I b — a 



(iv) Let T = R, (? : [0,1] — )• (0,oo) and h{t) = 1. Applying Theorem 11151 with the convex 
and continuous function / = — In on (0, 00), a = and 6 = 1, we get: 

In / g{s)ds > / ln{g{s))ds. 
Jo Jo 

Therefore, 



/ gis)ds > exp f / \n{g{s))dsj . 



(v) Let T = Z and n € N. Fix a = 1, b = n-\-l and consider a function g : {1, . . . , n + 1} — ?■ 
(0,cx3). Obviously, / = — In is convex and continuous on (0, 00), so we may apply 
Jensen's inequality to obtain 
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and hence 



In 



^ / n n+1 



t=i 

1 rn+l 

> - / ln{g{t))<)at 

n 



t=2 



In 



n 



n+l 



9{t)0at 



n+l 



t=2 



t=i t= 

' n "I ~ (n+l 

ln\llgit)\ +ln\llg{t) 



.t=i 



, i=2 



I / n n+l \ ( n -\- (n+l 

-{aY,9{t) + {l-a)Y,9m > \\{9m \\{9{t) 



t=i 



t=2 



a=i 



. t=2 



When a = 1, we obtain the well-known arithmetic-mean geometric-mean inequality: 



-, n ^ n 



(6.39) 



t=i 



.t=i 



When a = 0, we also have 



, n+l rn+l 



t=2 



. t=2 



(vi) Let T = 2^0 and iV e N. We can apply Theorem [IT5] with a = 1, 6 = 2^ and 5 : {2^ : 
0<k<N} ^ {0,00). Then, we get: 



In 



A 9{t)0at 

2^-1 



/ if5(t)At jfgmtX 

In < a \^r : \- [l — a) ^ 



2^-1 



2^-1 j 



1^ ^ «ELy2^g(2^) ^ (L- «) Eti 2^=5(2'=) 



> 



2^-1 

/f ln(g(t))0„t 

OA' _ 1 



2^-1 
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" 2^-1 + (^ - "^ 2^-1 

2^-1 2^-1 

2^-1 ^ 2^-1 



2^-1 2^-1 

I fc=o J U=i J 



fc=0 J Kk=l 

We conclude that 

1^J^ES^2^-) + (1 - a) Ef=i 2^^5(2'=) 1 



2^-1 



1 



fc=0 J Kk=l J 

On the other hand, 

N-l N N~l 

a Y, 2'5(2'=) + (1 - «) E 2'5(2'=) = ^ 2^g(2'=) + ag{l) + (1 - a)2^<7(2^). 

fc=0 fc=l fc=l 

It follows that 

Ef=i' 2^=5(2^) + ag{l) + (1 - a)2^5(2^) 



2^-1 






I fc=0 J K.k=l 

In the particular case when a = 1 we have 



2^-1 
and when a = we get the inequality 



fc=0 



11(5(2'))' 



%^.{n(.(.'r}'"- 
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6.3.1 Related Oa-integral inequalities 

1 1 
The usual proof of Holder's inequality use the basic Young inequality xpyi < - + ^ for 

nonnegative x and y. Here we present a proof based on the application of Jensen's inequality 

(cf. Theorem II 15p . 

Theorem 120 (Holder's inequality). Let T be a time scale, a, 6 € T with a < b, and f , g, 

h S C{[a,b]j, [0, oo)) with J h{x)g''{x)()a-c > 0, where q is the Holder conjugate number of 
p, i.e. ^ + ^ = 1 with 1 < p. Then, we have: 

h{x)f{x)g{x)0ax < (j h{x)fP{x)Oaxy (j h{x)g'i {x)0o.^' ■ (6.40) 

Proof. Choosing f{x) = x^ in Theorem 11151 which for p > 1 is obviously a convex function 
on [0, oo), we have 

( la\h{sMs)Oas Y ^ fjh{s)\{g{s)r0^s 

V Ia\Hs)\OaS J - J^\h{s)\OaS 

Inequality ()6.40p is trivially true in the case when g is identically zero. We consider two cases: 
(i) g{x) > for all x G [a, 6]t; (ii) there exists at least one t G [a, b]j such that g{x) = 0. We 

— q 

begin with situation (i). Replacing g hj fg p and |/i(x)j by hg'^ in inequality ()6.4ip . we get: 

l a Kx)a''{x)f{x)g p (x)OqX \ ^ /g Hx)9''{x){f{x)g p {x))P<)aX 
fl^h{x)gi{x)0ax ) ~ JaKx)gi{x)(>aX 

Using the fact that - + - = 1, we obtain that 

h{x)f{x)g{x)0ax < (j h{x)fP{x)0^xy (j h{x)gi{x)0^xy . (6.42) 

We now consider situation (ii). Let G = {x G [a, &]t | g{x) = 0}. Then, 

h{x)f{x)g{x)OaX = / h{x)f{x)g{x)OaX+ / h{x)f{x)g{x)OaX 



a,b]i-G JG 

h{x)f{x)g{x)()aX 



[a,h]^-G 

because /^ h{x) f {x)g{x)()aX = 0. For the set [a, b]j — G we are in case (i), i.e. g{x) > 0, and 
it follows from (|6.42p that 

h{x)f{x)g{x)OaX = / h{x)f{x)g{x)OaX 



\a,b]T-G 

\ - - 

<{ I h{x)fP{x)0^xf If h{x)g^x)(}axy 

la,b]r-G J \J[a,b]j-G J 

< ( / h{x)fP{x)Oax] ' ( [ h{x)g^x)^^x^ " 
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Remark 121. In the particular case h = 1, Theorem 11201 gives the Oa version of classical 
Holder's inequality: 



\fix)g{x)\OaX< ( [ \f\^ix)OaX 



\9\''ix)0aX] , 



where p > 1 and q = -^ . 

Remark 122. In the special case p = q = 2, (|6.40|) reduces to the following Oq Cauchy-Schwarz 
integral inequality on time scales: 



\f{x)g{x)\0ax < 



P{x)0^x] ( I g^{x)0ax 



We are now in position to prove a Minkowski inequality using Holder's inequality ()6.40p . 

Theorem 123 (Minkowski's inequality). Let T be a time scale, a, h ^ T with a < b, and 
p > 1. For continuous functions f,g : [a, b]j —^M we have 



l(/ + 5)(^)rOa.xV<f / |/(x)rOaX) +(/ \g{x)\POaX 



Proof. We have, by the triangle inequality, that 



\f{x)+g{x)\PO^x= / \f{x)+g{x)r'\f{x)+g{x)\0^x 

J a 

< I \f{x)\\f{x)+g{x)r'0ax+ I \gix)\\fix)+g{x)r'0ax. (6.43) 



Applying now Holder's inequality with q = p/{p — 1) to (|6.43p . we obtain: 



f\f{ 

J a 



x)+g{x)\P()c,x< 



+ 



\f{xWOaX 



\g{x)TOaX 



|/(x) + <7(x)|(P-l)'^0aX 



\f{x)+g{x)\^-^^'^O^x 



['ifixWOaX 
J a 



+ 



\gix)\POaX 



\f{x)+gix)\POo.x 



Dividing both sides of the last inequality by 

rb 



\f{x)+gix)\POaX 



we get the desired conclusion. 



D 
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6.4 Some integral inequalities with two independent variables 

We establish some linear and nonlinear integral inequalities of Gronwall-Bellman-Bihari type 
for functions with two independent variables. Some particular time scales are considered as 
examples and we use one of our results to estimate solutions of a partial delta dynamic 
equation (cf. ()6.62p below). 

As was mentioned earlier, several integral inequalities of Gronwall-Bellman-Bihari type on 
time scales but for functions of a single variable were obtained in the papers [2], [8], [83] and 
|93j . However, to the best of the author's knowledge, no such a paper exists on the literature 
when functions of two independent variables are considered. Therefore we aim to give a first 
insight on such a type of inequalities. 

6.4.1 Linear inequalities 

Throughout we let Ti = [ai, oo)ti and T2 = [a2, oo)t2) for ai € Ti, 02 S T2 being Ti and T2 
time scales. 

Theorem 124 (cf. 05]). Let u{ti,t2),a{ti,t2), f{ti,t2) G C(ti x t2,]R^) with a{ti,t2) 
nondecreasing in each of the variables. If 

U{ti,t2)<a{ti,t2)+ / f{si,S2)u{si,S2)AiSiA2S2, (6.44) 

J ai J a2 

for{ti,t2) € ti X ¥2, then 

u{ti,t2) < a{ti,t2)ejt2j^^^,.^)^^^^{ti,ai), (ti,t2) S ti x ¥2. (6.45) 

Proof. Since 0(^1,^2) is nondecreasing on (^1,^2) € Ti x T2, inequality ()6.44p implies that, for 
an arbitrary e > 0, 

rti rt2 

r{ti,t2)<l+ / f{si,S2)r{si,S2)AiSiA2S2, 

where r{ti^t2) = ^a Vvfe * Define ^(^1,^2) by the right-hand side of the last inequahty. Then 

^ (^^Xj^) = /(*i'*2)r(ti,t2) < f{tiM)v{ti,t2), {tiM) G tj X t^. (6.46) 

From ()6.46p and taking into account that f(ti,t2) is positive and nondecreasing, we obtain 



V{tl,t2)v{ti,a2{t2)) 

from which it follows that 



<f{tiM), 



\ 1, ^JAiti \ Aiti J ^ j,^^ ^ ^ , Aiti A2t2 



<f{tl,t2) + 



V{ti,t2)v{ti,a2{t2)) ' V{ti,t2)v(ti,a2{t2))' 
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The previous inequality can be rewritten as 

n / dv{tl,t2) 

" { Alii 



A2t2 \v{tut2) 



<f{tl,t2)- 



dv{ti,t2) 



Integrating with respect to the second variable from 02 to ^2 and noting that ^^^' ^' |(ti,a2)~ 0; 
we have 



that is, 



v{t 



^-^< r f{tl, 82)^282, 
1)I2J Ja2 



%^< r f{tl,S2)A2S2vih,t2). 



Fixing ^2 £ Tr2 arbitrarily, we have that p{ti) '■= f ^ f{tii -52)^252 G '^''' and by Theorem 



0.2 



v{ti,t2) < ep{ti,ai), 
since v{ai,t2) = 1. Now, inequality (j6.45p follows from the inequality 

u{ti,t2) < [a{ti,t2) + e]v{ti,t2), 
and the arbitrariness of e. 



D 



Corollary 125 (cf. Lemma 2.1 of ^65j). Let Ti = T2 = M and assume that the functions 
u{x,y),a{x,y), f{x,y) € C([xo,oo) x [yo,oo),IR^) with a{x,y) nondecreasing in its variables. 

If 

(6.47) 



u{x,y) < a{x,y) + / / f{t,s)u{t,s)dtds 

Jxii Jyo 

for {x,y) G [xo,oo) x [yQ,oo), then 



u{x,y) <a{x,y)exp[ I I f{t,s)dtds], 

'xo Jyo 



X py 



(6.48) 



for {x,y) e [xo,oo) x [yo,oo). 

Corollary 126 (cf. Theorem 2.1 of [88]). Let Ti = T2 = Z and assume that the functions 
u{m, n), a{m, n), f{m, n) are nonnegative and that a{m, n) is nondecreasing for m € [moi oo)n 
Z and n G [no, 00) n Z, niQ^UQ € Z. // 

m—l n—1 

u{m,n) < a{m,n) + N^ N^ /(s, f)n(s,i), (6.49) 



=mo t=no 



for all (m, n) € [niQ, 00) n Z x [no, 00) PI Z, then 



m—l 



u{m,n) < a{m,n) I I 



=mo 



1 + E /("'*) 

t=no 



.50) 



for all (m, n) € [mo, 00) n Z x [no, co) PI Z. 
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Remark 127. We note that, following the same steps of the proof of Theorem 11241 it can be 
obtained other bound on the function u, namely 

u{ti,t2) < a{ti,t2)ejHf^,^t2-^^^,^{t2,a2). (6.51) 

When Ti = T2 = M, then the bounds in ()6.5ip and ()6.45p coincide (see Corollary I125p . But if, 
for example, we let Ti = T2 = Z, the bounds obtained can be different; moreover, at different 
points one can be sharper than the other and vice-versa. Note up the following example: 

Let /(ti,t2) be a function defined by /(0,0) = 1/4, /(1,0) = 1/5, /(2,0) = 1, /(0, 1) = 
1/2, /(1, 1) = and /(2, 1) = 5. Set ai = 02 = 0. 

Then, from (j6.45p we get 

n(2,l)<a(2,l)^. 



while from ()6.5ip we get 



147 
n(3,2)<a(3,2) — , 



29 

n(2,l)<a(2,l)-, 

n(3,2)<a(3,2) — . 



We now present the particular case of Theorem 11241 when Ti = Z and T2 = M. 
Corollary 128 (cf. [IS]). Let T\ = TL and T2 = M. Assume that the functions u{t,x), a{t,x) 



and f(t,x) satisfy the hypothesis of Theorem, ] 124\ for all {t,x) € Ti x T2, with ai = 02 = 0. // 

u{t,x) <a{t,x) + y2 f{k,T)u{k,T)dT, (6.52) 

for all {t,x) G Ti X T2, then 



^ V fx 

u{t,x) <a{t,x)Y{ 1+/ fik,T)dT 

I — n L JQ 



k=0 

for all{t,x) e ti X ¥2. 



(6.53) 



We now generalize Theorem 11241 If in Theorem 11291 / does not depend on the first two 
variables, then we obtain Theorem 11241 

Theorem 129 (cf. [IS]). Letu{ti,t2),a(ti,t2) S C(Tri x ']r2,M(|'), with a{ti,t2) nondecreasing 
in each of the variables and f{ti,t2, 81,82) € C(S', MjJ') be nondecreasing in ti and t2, where 
S = {{ti,t2,8i,82) e ti x t2 X ti X ¥2 :ai <si <ti,a2 < S2 < t2}. // 

rti rti 

U{ti,t2) < a{ti,t2) + / f{ti,t2,8i,82)u{8i,82)Ai8iA282, (6.54) 

J ai J a2 
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for{ti,t2) € ti xfs, then 

u{ti,t2) < a(*i'*2)ej^*2 j(j^^j2,ii,s2)A2S2(*i'"i)' (*i'*2) G ti X ¥2. (6.55) 

Proof. We start by fixing arbitrary numbers tj G Ti and ^2 ^ ^'2 and consider the function 
defined on [ai, t*] n ti x [02, ^2] n ¥2 by 

v{ti,t2) = a{tl,t2) + e+ / /(ti,t2,si,S2)^i(si,S2)AiSiA2S2, (6.56) 

J ai J a2 

for an arbitrary e > 0. From our hypothesis we see that 

u{tut2)<v{h,t2), for all (ti,t2) G [ai,tt]nti X [02,^2] nt2. (6.57) 

Moreover, A-differentiating with respect to the first variable and then with respect to the 
second, we obtain 



A2i2 V Alii 

< f{tl,t*2,ti,t2)v{h,t2), 

for all (ti,t2) G ['^^ij^il'^ n Ti X [02,42]" n T2. From this last inequality, we can write 






< /(*1;*2)*1)*2), 



V{ti,t2)v{ti,a2{t2)) 

hence, 

nn,t2)A2t2 [ Aiti ; ^Y7^ — s^t^ 



<f{t*l,t*2,ti,t2) + 



v{h,t2Hti,a2{t2)) -•'^i'2'1'^^ w(tl,t2)^^(tl,fT2(i2))' 

The previous inequality can be rewritten as 

dv{ti,t2) \^ 

< /(il)i2)*l)*2)- 



d I Aiti 



A2t2 \v{ti,t2) 



dvjtiM) 



A-integrating with respect to the second variable from 02 to ^2 and noting that ^^^' '^' \[ti,a2)- 
0, we have 

-t^.< / f{t\AM.S2)/^2S2. 
vitiM) Ja2 

that is, 

5^;(ti,t2) / Z"*^ 



Aiti 



< f' fitl,t*2,h,S2)A2S2viti,t2). 



'a2 

Fix t2 = ^2 and put p(ti) := /J /(t^ t^, ti, S2)A2S2 G 7^+. By Theoremge] 

t'(ii,i2) < (a(ti,t2) + e)ep(ti,ai). 

Letting ti = tj in the above inequality, and remembering that t^, ig and e are arbitrary, 
inequality (|6.55p follows. D 
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6.4.2 Nonlinear inequalities 

Theorem 130 (cf. 05]). Let u(ti, t2), /(^i, ^2) G C(ti x t2,]R^). Moreover, let 0(^1,^2) G 
C(Ti X T2,M'^) and be a nondecreasing function in each of the variables. If p and q are two 
positive real numbers such that p > q and if 

U^{tl,t2)<a{ti,t2)+ r /"'/(si,S2K(si,S2)AiSiA2S2, (6.58) 

for{ti,t2) G ti xf2, then 



u{ti,t2) < ap{ti,t2) 



e „, £-1 (ti,ai) 



, (ti,i2)Gtixt2. (6.59) 



Proof. Since 0(^1,^2) is positive and nondecreasing on (^1,^2) G Ti x T2, inequality (|6.58p 
implies that, 

uP{h,t2)<a{ti,t2)(l+ r r f{si,S2)^^-j^^AiSiA2S2). 
V Jai Ja2 a{si,S2) J 

Define v(ti,t2) on Ti x T2 by 

V{t,,t2) = l+ r r f{si,S2)^^-j^^A^S,A2S2. 
Ja^ Ja2 a{Si,S2) 

Then, 

d fdv{ti,t2)\ ,,, , .u^^itiM) ^ f,. . . a.^1 . . 3. 

IT— — r-7 = f{ti,t2)^-—— < f{tl,t2)av {ti,t2)vP ti, ta • 

A2t2 V Aiti ) a{ti,t2) 

q_ 

Now, note that vp{ti^t2) < ^(^1,^2) ^^d therefore 

^(5^)<m,«.)»5-'(ti,feM«i.<.). 

From here, we can follow the same procedure as in the proof of Theorem 11241 to obtain 

v{ti,t2) < ep{ti,ai), 

where p{ti) = JJ^ /(ti,S2)«^~ {ii,S2)A2S2. Noting that 

1 1 

u{ti,t2) < ap{ti,t2)vp{ti,t2), 

we obtain the desired inequality (|6.59p . D 

Theorem 131 (cf. 08]). Letu{ti,t2),a{ti,t2) G C{TixT2,M.q), witha{ti,t2) nondecreasing 
in each of the variables and f{ti,t2, si, S2) G C(5, M^), where S = {(ti, ^2, si, S2) G Ti x T2 x 
Ti X T2 : fli < si <ti,a2<:S2< ^2}- If P '^^^ Q '^^s two positive real numbers such that p> q 
and if 

uP{ti,t2)<a{ti,t2)+ f f /(ii,t2,si,S2)n«(si,S2)AiSiA2S2, (6.60) 
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for all (^1,^2) € Ti x T2, then 

1 
u{ti,t2) < ap{ti,t2) 



2 „, 9_i (ti,ai) 

Jl^aV \ti,S2)f{tl,t2,tl,S2)A2S2^ " " 

for all {h,t2) € ti X ¥2. 

Proof. Since a{ti,t2) is positive and nondecreasing on (ti,t2) G Ti x T2, inequality ()6.60p 
implies that, 



u''{h,t2)<a{ti,t2)(l+ r f' f{tl,t2,Si, 



U''{SI,S2) . . , 

S2)^ -A1S1A2S2 I • 

a{si,S2) 



Fix t* € Ti and tg S ^2 arbitrarily and define a function v{ti,t2) on [ai, t*] fl Ti x [02, ^2] ^ '^2 

by 

U'^{SI,S2) 



V{t,,t2) = l+ r [' f{tl,t*2,SuS2)^^-^^^A,SiA2S2. 
Jai Ja2 a{Sl,S2) 



Then, 



d fdv{ti,t2) 



A2t2 



f dv{ti,t2) \ 



f{tl,t2,tl,t2) 



a{si,S2) 



a{ti,t2) 
< f{t\, t*2, h,t2)ap-\ti,t2)vv (ti, ta). 



Since t'p(ti,t2) < ^(^1,^2)? we have that 

^ (^^1^) < f{tl,t*2,h,t2)J~\t,,t2)vih,t2). 

We can follow the same steps as before to reach the inequality 

^^^^1^ < r f{tl,t*)g{tl,t*2,tl,S2)av~\h,S2)A2S2v{h,t2). 

Fix t2 = t*2 and put p{h) := f^l f{tl,q)g{tl,tlti,S2)ap~\tu 82)^282 G 7^+. Again, an 
application of Theorem UH] gives 



f(ti,t2) < ep{ti,ai), 
and putting ti = tl we obtain the desired inequality. 



D 



We consider now the following time scale: let {afc}fcgN be a sequence of positive numbers 
and let t^ £ R. Let 



ifc=io+E""' ^^^' 



n=l 



and assume that limfc^oo tk = co- Then we can define the following time scale T" = {f^ : k G 
No}. Further, for p € 7^ (cf. [H Example 4.6]), 

A,- 
ep{tt,to) = n (^ + o^nP{tn-i)), for ah k G Nq. (6.61) 



n=l 
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Now, for two sequences {ofc, /3fc}fceN and two numbers tg , tg € M as above, we define two time 
scales T" = {t^ : k S No} and T'^ = {t^ : k G No}- We can therefore state the following 
result: 

Corollary 132 (cf. [E]). Let u{t,s),a{t,s) be defined on T" x T° , be nonnegative with a 
nondecreasing in its variables. Further, let f(t,s,T,^), where (t,s,T,^) € T" x T^ x T" x T^ 
with T < t and S, < s, be a nonnegative function and nondecreasing in its first two variables. 
If p and q are two positive real numbers such that p > q and if 



w 



'{t,s)<a{t,s)+ Yl E /^"(T)/(0/(i,^,T,eK(r,0, 



for all {t,s) G T" X T'^, where 11°' and jj,^ are the graininess functions o/T" and T^ , respec- 
tively, then 



u{t,s) < ap{t,s) 
for all {t, s) G T" X T'^, where e is given by 116. 6 1\) 



^f;pf{t,s)ap-\t,Og{t,s,WAP^^ ' 0' 



We end this section showing how Theorem 11301 can be applied to estimate the solutions 
of the following partial delta dynamic equation: 



d (du^{ti,t2) 



F{ti,t2,uitut2)), 



A2i2 V Alii 

with the given initial boundary conditions (we assume that ai = 02 = 0) 

u\ti,0) = g{ti), n2(0, t2) = h{t2), g{0) = 0, /i(0) = 0, 



(6.62) 



(6.63) 



where F G C(ti x t2 x M(]",M(]"), g G C(ti,M+), h G C(t2,M(J'), with g and h nondecreasing 
functions and positive on their domains except at ai and 02, respectively. 

Theorem 133. Assume that on its domain, F satisfies 

F{ti,t2,u) < t2U. 

If u{ti,t2) is a solution of the IBVP above for (^1,^2) G Ti x T2, then 



u{tiM)< Vi9{ti) + h{t2)) 



''lo^S2ig{h)+h{s2)rhA2S2^^^'^^ 



(6.64) 



for (^1,^2) G Ti x T2, except at the point (0,0). 



Proof. Let u{ti,t2) be a solution of the given IBVP (|6.62|) - (j6.63p . Then it satisfies the 
following delta integral equation 

rti rt2 

U^{tl,t2) = g{ti)+h{t2)+ / F{si,S2,u{si,S2))AiSiA2S2. 

Jo Jo 
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The hypothesis on F imply that 

rti f-t2 

U^itl,t2) < g(ti) + h{t2)+ / S2u{si,S2)AiSiA2S2. 

Jo Jo 
An appHcation of Theorem ll30l with a(ti . U) = g{ti) + h{t2) and f{ti,t2) = t? gives (|6.64p . D 

6.5 State of the Art 

The results of this chapter are published or accepted for publication in [Til SSI HSl HE]. 
Moreover the three papers |4H H2| H7] are already available in the literature and contain 
related results. Some partial results of [H] were presented in the Time Scales seminar at the 
Missouri University of Science and Technology. 

The study of integral inequalities on time scales (with both A-integrals or Oo-integrals) 
and its applications is under strong current research. We provide here some more (recent) 
references within this subject: [T2l [T3| [26| [361 [371 [68] . 
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Chapter 7 



Discrete Fractional Calculus 



We introduce a discrete-time fractional calculus of variations on the time scale /iZ, h > 0. 
First and second order necessary optimality conditions are established. Examples illustrating 
the use of the new Euler-Lagrange and Legendre type conditions are given. They show that 
solutions to the considered fractional problems become the classical discrete-time solutions 
when the fractional order of the discrete-derivatives are integer values, and that they converge 
to the fractional continuous-time solutions when h tends to zero. Our Legendre type condition 
is useful to eliminate false candidates identified via the Euler-Lagrange fractional equation. 

7.1 Introduction 

The Fractional Calculus is an important research field in several different areas: physics 
(including classical and quantum mechanics as well as thermodynamics), chemistry, biology, 
economics, and control theory [77t I91| . It has its origin more than 300 years ago when 
L'Hopital asked Leibniz what should be the meaning of a derivative of non-integer order. 
After that episode several more famous mathematicians contributed to the development of 
Fractional Calculus: Abel, Fourier, Liouville, Riemann, Riesz, just to mention a few names. 
In [76] Miller and Ross define a fractional sum of order z^ > via the solution of a linear 
difference equation. They introduce it as 

^-■^fit) = ^ E(t - ct(.))(^-^)/(s). (7.1) 

^ ' s=a 

Definition (j7.ip is the discrete analogue of the Riemann-Liouville fractional integral 



„D-V(x) = ^ J\x - sr-^fis)ds 



of order v > 0, which can be obtained via the solution of a linear differential equation |761 l77j. 
Basic properties of the operator A~'^ in ()7.ip were obtained in [76]. More recently, Atici and 
Eloe introduced the fractional difference of order a > by A"/(t) = A(A~(^~")/(t)), and 
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developed some of its properties that allow to obtain solutions of certain fractional difference 
equations [151 US] • 

The fractional differential calculus has been widely developed in the past few decades due 
mainly to its demonstrated applications in various fields of science and engineering (BSJ [771 
I86j . The study of necessary optimality conditions for fractional problems of the calculus of 
variations and optimal control is a fairly recent issue attracting an increasing attention - see 
in El [Ml [Ml [221 [S3] and references therein - but available results address only the continuous- 
time case. Therefore, it is pertinent to start a fractional discrete-time theory of the calculus 
of variations, namely, using the time scale T = hZ, h > 0. Computer simulations show that 
this time scale is particularly interesting because when h tends to zero one recovers previous 
fractional continuous-time results. 

Our objective is two- fold. On one hand we proceed to develop the theory of fractional 
difference calculus, e.g., introducing the concept of left and right fractional sum/difference (cf. 
Definitions 1 1 38 1 and 1 1 44 p . On the other hand, we believe that the results herein presented will 
potentiate research not only in the fractional calculus of variations but also in solving frac- 
tional difference equations, specifically, fractional equations in which left and right fractional 
differences appear [see the Euler-Lagrange Equation (|7.15p ]. 

Main results so far obtained by us contain a fractional formula of /i-summation by parts 
(Theorem I146p , and necessary optimality conditions of first and second order (Theorems 11511 
and 11541 respectively) for the proposed /i-fractional problem of the calculus of variations 
(j7.14p . We also provide some illustrative examples. 

We remark that it remains a challenge how to generalize the present results to an arbitrary 
time scale T. This is a difficult open problem since an explicit formula for the Cauchy function 
(cf. Definition [27]) of the linear dynamic equation y = on an arbitrary time scale is not 
known. 



7.2 Preliminaries 

It is known that y{t, s) := Hn-i{t, (t(s)) is the Cauchy function (cf. Definition I27p for y = 
[241 Example 5.115], where -ffn-i is the time scale generalized polynomial /i„_i of Definition 
[TSl (we use here H instead of h in order to avoid confuse with the parameter h of the time 
scale T = /iZ). 

From now on we restrict ourselves to the time scale T = /iZ, h > 0. 

If we have a function / of two variables, f{t,s), its partial forward /i-differences will 
be denoted by At^h and Ag^h, respectively. We will make use of the standard conventions 
EtZl m = 0, c G Z, and U^, f{i) = 1. 

Before giving an explicit formula for the generalized polynomials H^ on T = KL we 
introduce the following definition: 
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Definition 134 (cf. [E]). For arbitrary x,y G M the /i-factorial function is defined by 

Jy) ._ uy r(f + l) 

' ■" r(| + i-y)' 

where T is the weh-known Euler gamma function, and we use the convention that division at 
a pole yields zero. 

Remark 135. For h = 1, and in accordance with the previous literature [see formula (j7.ip ]. 
we write x^^' to denote xf^ . 

Proposition 136. For the time-scale T one has 

it - si^^ 

HJt, s) := — ^^ for all s,t G hZ and /c G Nq . (7.2) 

k! 

To prove ()7.2p we use the following technical lemma. We remind the reader the basic 

property r(x + 1) = xT{x) of the gamma function. 

Lemma 137. Let s E T. Then, for all t £T'^ one has 

f (t-.r" l ('-»)!f' 

'■'• \ (fc + 1)! j " ki. ■ 

Proof The equality follows by direct computations: 

( (t - st-^'^ 1 _ 1 / (ait) - s)^'^ it - sfr^ \ 



A 



t,h 



{k + i)\ j h] (k + iy. {k + iy. j 

/I'^+i r T{{t + h-s)/h + l) T{{t-s)/h + l) 



h{k + 1)! i r((i + h- s)/h + 1 - (fc + 1)) r((t - s)/h + 1 - (A; + 1)) 
h^ r T{{t-s)/h + 2) T{{t-s)/h + l) 



{k + l)\ \T{{t-s)/h + l-k) T{{t-s)/h-k) 

h^ \{{t-s)/h + l)T{{t-s)/h + l) V{{t-s)/h + l) 



{k + 1)! i ((t - s)/h - k)T{{t - s)/h - k) T{{t - s)/h - k) 

h^ r {k + l)T{{t-s)/h + l) \ _h'' { T{{t-s)/h + l) 



{k + 1)! i {{t - s)/h - k)T{{t - s)/h - k) ] k\ \ T{{t - s)/h + l-k) 
k\ 



Proof, (of Proposition I136p We proceed by mathematical induction. For k = 



Assume that ()7.2p holds for k replaced by m. Then, by Lemma 11371 



ft ft I \V'"'/ 

H.m+i{t,s) = / H^{t,s)/\t= / ^^"7^ At 



D 



(m + 1)! 
which is ()7.2p with A: replaced by m + 1. D 
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Let 2/1 (t), ... , Unit) be n linearly independent solutions of the linear homogeneous dynamic 
equation y = 0. From Theorem 1281 we know that the solution of ()1.12p (with L = A" and 
to = a) is 

y{t) = A~-m= / ^' yl^ /(.)A. = -i- J2{t-a{kh)t-''^f{kh)h. 

'''' ^ ^"^^ ^ ^'^> k=a/h 

Since y "-{a) =0, i = 0, . . . ,n — 1, then we can write 

t/h—n 

A-"/W = ^ E {t-c7{kh)t-'^f{kh)h 

^"^^ k=a/h (7.3) 

Note that function t — ?• (A~"/)(i) is defined for t = a + n/i mod(/i) while function t ^ f{t) 
is defined for t = a mod(/i). Extending ()7.3p to any positive real value v, and having as an 
analogy the continuous left and right fractional derivatives [77], we define the left fractional 
/i-sum and the right fractional /i-sum as follows. We denote by J-'j the set of all real valued 
functions defined on the time scale T. 

Definition 138. Let / € J^j. The left and right fractional /i-sum of order i^ > are, 
respectively, the operators a^^'^ '■ J^T — ^ -^f + and h'^b'^ ■ ^T — > -^f- ■, Tj^ = {a ± vh : a <^ T}, 
defined by 

--U 

a^l'fit) = ^ £^* ''\t - a(.))i'^-^V(^)As = J-^ ^jt - <y{kh)t-'^f{kh)h 

h 
b_ 

Lemma 139. Let v > Q he an arbitrary positive real number. For any t G T we have: (i) 
lim,^OaA^^/(t + i./i) = f{t); (ii) lim^^oh^-'^fit - uh) = f{t). 

Proof. Since 

1 Mt) 



1 r^^' I 1^ 

„A-^/(t + z./i) = ^ y^ (t + ^/i - a(s))t'^-^V(s)A. 



h 

TM 



h 






^^*^ + ThTTT) E (^ + ^^ - <^(.kh)fr'^f{kh)h , 



k=t 
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it follows that lim^^oaA^j'^/(i + i^h) = f{t). To prove (ii) we use a similar method: 



^'^(*^ + foTTTT ^ (^/^ + ^/^-^W)r'V(fc/^)/i 



T{u + 1) 



."(t) 



and therefore lim^^o h'^i, "^ fit — vh) = f{t). 

For any t € T and for any z^ > we define aA^/(t) := /iA|J/(t) := /(t) and write 
aA-'^f{t + Z./1) = h'^fit) + ^.^-^ y^ (t + i^h- a{s))t'^f{s)As , 

nAl^fit) = h^f{t - uh) + —-— / {s + uh- a{t))t'^fis)As . 
Theorem 140. Let f G J^j and u >0. For all t eT'^ we have 



D 



(7.4) 



^A-''f^{t + uh) = {aA-''f{t + uh)y 



.{t + uh-at-'^fia). (7.5) 



r(z. + l) 
To prove Theorem 11401 we make use of the following 

Lemma 141. Let t G T". The following equality holds for all s G T".- 

A,,^,(^it + uh-s)i'-'^f{s))) 

= {t + vh- ais))i'-'^f%) -{v- l)(t + vh- a{s)ir^^f{s) . (7.6) 
Proof. Direct calculations give the intended result: 
A,^f,({t + uh-si;^~^'^f{.s, 

)h 

it + ^h-a{s))t-''^-{t + yh-st-^'^ .,.^,^^, , .. (.-1) 
fKs)^{t^vh-a{s))\ 



As,h ({t + uh- sf;:"" ) f{s) + {t + vh- a{s))t-''^ f^{s) 



'U~l 



h 

-p f t+vh-a(s) 1 



^ ( t+uh-a(s) ^ ^ _ ^^ _ ^^ 



m 

h 

+ {t + uh-a{s)fr'^ f'^is) 

fis 



'U-1 



p I t+vh—s 



+ 1) 



■t+uh- 



+ l-(^-l)) 



h 

f{sW~' 
f{sW-^ 



v-2 



Y-( i-\-vh—s \ T^t t+vh—s \ -\\ 



r(i^ + i) r(^ + 2) 



+ (t + ^/,_^(5))(^--l)/A(^ 



r(^±x^) (t 
r(¥ + 2) 



h 



+ 1 



t + i^h — s 
h 



+ it + uh-a{s))t'^f^{s) 



■pf t+uh—s—h _i_ ■\\ 



,(1.-2) 



{iy~l) fA, 



(j, _ i)(i + ^/, _ a(s))r 7(«) + it + '^h- a{s))\r" f^{s) 
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where the first equality follows directly from (|1.3 



D 



Remark 142. Given an arbitrary t S T'^ it is easy to prove, in a similar way as in the proof 
of Lemma I14H the following equality analogous to ()7.6p : for all s € T*^ 

= (^ _ l)(s + ^/i - a{t))l:'-'^r{s) + is + ,.h- a(i))r'V^(s) . (7.7) 
Proof, (of Theorem I140p From Lemma 11411 we obtain that 

aA-'^f'^it + uh) = h'^f^it) + ^;^^ j\t + vh- a{s))t'^f'^{s)As 



h'f'^it) + 



T{u + 1) 



{t + uh-s)t'^f{s) 



r(i^ + i)A 



<x{t) 



{u-2) 



{i^-m + i^h-a{s))):;-'>fis)As 



(7.8) 



T(u + l) 



-f{a) + h-f'^{t) + vh''-^f{t) 



+ 



r(i^ + i)7a 

We now show that {aAj^^fit + uh))^ equals dLH]): 



{u-l)it + i^h-ais)fr^^f{s)As. 



{,A-'^f{t + uh))^ = - 



-h'fit) 



h-f{a{t)) + 



n^+^)ja 



a(t) 



{ait)+,.h-a{s))^r'^f{s)As 



.I'fA 



h'r{t) + 



T{v + l) 



hViv + 1) 



\t + vh-a{s)t^'^f{s)As 



{u-l) 



{<y{t) + vh-a{s)f;:-'>f{s)As 



+ h^~'uf{t) 



J a 

fHt) + Y^^^ £ A,,, ((i + uh- C7(.))i^-^)) f{s)As + h'^~'uf{t) 



D 



Follows the counterpart of Theorem 11401 for the right fractional /i-sum: 
Theorem 143. Let f G Tj and u > 0. For all t G T" we have 



.A-^.^r (t - uh) = w^;^Yy(b + ^h- a{t))r" f{h) + (.A-^/(t - vh)Y 



(7.9) 
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Proof. From ()7.7p we obtain that 

,A-^)/^(t - uh) = h-^f^it) + — ^ /' {s + uh- a{t))t'^fHs)As 

^ \^ ' ^) J a(t) 



h^rit) + 



uis + uh-a{t))^;^~^^ 



fis) 



s=b 



S = CT{t) 






(7.10) 



We show that {hA~''f{t - vh))^ equals (17111]) : 

/■o"(^) 



h''r{t) + 



/ir(i/ + 1) 






vh--'f{a{t)) 



- / (s + z./i-a(t))t'^-^V(5)As 

r(i^ + i) y<x2(t) ^ ^ 



h^rit) 



u 



^{^+l)Ja 



(t) 
a(b) 



2(.\\{''-'^) 



■V j-A 



h-rit) 



r(i^ + i)y. 



2(t) 



(^ - 1)(5 + ^/i - a\t)X- l^s)As - v\f~'!{a{t)) 



Av-2) 



it) 



(^ - l)(s + Z./1 - a(t))r 7(*)^^ - ^h-~'f{^{t)) 



U 



Definition 144. Let < a < 1 and set 7 := 1 — a. The left fractional difference aA^/(t) 
and the right fractional difference hA'^f{t) of order a of a function / € Fj are defined as 



,A^/(i):=(aA-^/(t + 7/i))^ and ,A^/(t) := -UA-^/(i - 7/^))' 



for all t € T''. 
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7.3 Main results 

Our aim is to introduce the /i-fractional problem of the calculus of variations and to prove 
corresponding necessary optimality conditions. In order to obtain an Euler-Lagrange type 
equation (cf. Theorem I15ip we first prove a fractional formula of /i-summation by parts. 



7.3.1 Fractional /i-summation by parts 

A big challenge was to discover a fractional /i-summation by parts formula within this time 
scale setting. Indeed, there is no clue of what such a formula should be. We found it eventually, 
making use of the following lemma. 

2 

Lemma 145. Let f and k be two functions defined on T'^ and T'^ , respectively, and g a 

2 

function defined on T" x T" . Then, the following equality holds: 



fit) 



g(t, s)k{s)As 



At 



p{b) 



kit) 



a{t) 



g{s,t)fis)As 



At. 



Proof. Consider the matrices R = [f{a + h),f{a + 2h), ■ ■ ■ , f{b — h)], 



Ci 



g{a + h,a)k{a) 
g[a + 2h, a)k{a) + g{a + 2/i, a + h)k[a + h) 

g{b - h, a)k{a) +g{b-h,a + h)k{a + h) -\ \- g{b - h,b - 2h)k{b - 2h) 



Co 



g{a + h, a) 
g{a + 2/i, a) 



C. 



g{b-h,a) 
Direct calculations show that 





g{a + 2/i, a + h) 

g{b — h,a + h) 



, C4 













_g{b-h,b- 


-2h) _ 



, r „J -, b/h-1 i-1 

fit) J g{t,s)k{s)As At = h^ Yl /(^^) Z 9{ih,jh)k{jh)=h^R-Ci 

i=a/h j=a/h 

= h^R ■ [k{a)C2 + k{a + h)C-i + ■ ■ ■ + k{b - 2h)C4 

b/h-l b/h-l 

= h^ k{a) J2 9{jh,a)f{jh)+k{a + h) J] g{jh,a + h)f{jh) 

j=a/h+l j=a/h+2 

b/h-l 

+ ... + k{b-2h) J2 9{jh,b-2h)f{jh) 
j=b/h~i 

b/h-2 b/h-l ^^(fe) 



Y, kiih)h Y 9ijh,ih)fijh)h= kit) / gis,t)f{s)A 

,, ■ /■! \ I, J a J ait) 



i=a/h 



j=a{ih)/h 



At. 
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D 

Theorem 146 (fractional /i-summation by parts). Let f and g be real valued functions defined 
on T'^ and T, respectively. Fix < a < 1 and put 7 := 1 — a. Then, 

/ f{t)aKg{t)^t = hy{p{h))g{h) - hy{a)g{a) + / ;,A°(,)/(t)ff'^(t)Af 

J a J a 

+ TX^^^''^ [ly + ^^ " ")^'''^^(*)^* - [^ i^ + lh- ^(a))i""'V(i)At j . (7.11) 
Proof. By ()7.5p we can write 



/(t),A^5(i)At = / /(i)(aA-^5(i + 7/i))^Ai 

J a 



fit) 



,A-V(t + 7M + 



7 



(7-1). 



r 



^^(t + 7/.-a)rMa) 



At 



(7.12) 



J' /(t),A- V(t + 7/^)At + I' f^:^^it + lh- a)^^-'^f{t)g{a)At. 



Using (|7.4p we get 

[ f{t)aA-'g^{t + jh)At 

J a 

= [ fit) [/iV(i) + Y^:^^ f\t + lh- a{s)i;^-'\^{s)As 

= W I' f{t)g^{t)At + ^.^-^ I' fit) j\t + ^h- a{s))^r'^g^{s)AsAt 

= h^ I' fit)g^{t)At + ^.^-^ 1'^'^ g^{t) 1'^ {s + jh- cT(t))i7-^V(^)A.At 



At 



Pib) 



hy{p(h))[g{h) - g{p{h))] + / ff^(t),A-y (t - 7/i)At, 



where the third equaUty follows by Lemma 11451 We proceed to develop the right-hand side 
of the last equality as follows: 

<-p{b) 

'^^' ^^P{b)' 

t=p{b) 



hy{pib))[g{b) - g{p{b))] + I g\t)HA~:<f{t - ^h)At 



t=a 



= hyipib))[gib)-gipib))] + ^<7(i)hA-y(t-7/i)_ 

rpib) 

- / 5'^(t)(,A-y(t-7/i))^At 

Ja 

= hy{p{b))g{b) - h?f{a)g{a) 

- Y^^f^^^ [^ is + lh- aia))i'-'^f{s)As + J' ' (.A°(,)/(t)) <7'^(t)At, 
where the first equality follows from ()1.7p . Putting this into ()7.12p we get ()7.1ip . 



D 
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Let us prove that Theorem 11461 indeed generahzes the usual summation by parts formula 
Corollary 147. Suppose that h = a = 1 in Theorem \146\ Then, the next formula holds: 

/ fit)Agit)At = f{p{h))g{b) - f{a)g{a) - / Af{t)g%t)At. (7.13) 

J a J a 

Remark 148. Since /i = 1, then (|7.13p becomes 

6-1 fe-2 

Y, f{t)[g{t + 1) - g{t)\ = fib - l)g{b) - f{a)g{a) - ^[/(t + 1) - f{t)]g{t + 1). 

t=a t=a 

If we allow / to be defined at t = 6, then we have 

6-1 b-2 

^ /(t)b(t + 1) - g{t)] = fib - l)g{b) - f{a)g{a) - ^ifit + 1) - f{tMt + 1) 

t=a t=a 

= f{b)g{b) - f{a)g{a) - [f{b) - f{b - 1Mb) 

fe-2 

-Y.[f{t + i)-f{tMt + i) 

t=a 

b-1 



f{b)g{b) - f{a)g{a) - ^[/(t + 1) - f{t)]g{t + 1). 



t=a 

This is the usual summation by parts formula. 

7.3.2 Necessary optimality conditions 

We begin to fix two arbitrary real numbers a and (3 such that a,/3 £ (0, 1]. Further, we put 

7 := 1 — a and z> := 1 — /?. 

Let a function L{t, u,v,w) : T*^ x M x M x M — )• M be given. We consider the problem of 

minimizing (or maximizing) a functional C : J-j — )■ M subject to given boundary conditions: 

l-b 
£[y{-)] = / L{t, y%t), aA'^y{t),hA^,y{t))At -^ min, y{a) = A, y{b) = B . (7.14) 

J a 

Our main aim is to derive necessary optimality conditions for problem ()7.14p . 
Definition 149. For / € J-r we define the norm 

11/11 = max ir (t)| + max \aAlf{t)\ + max UA^ /(t)|. 

A function y G Fj with y{a) = A and y{b) = B is called a local minimum for problem (j7.14p 
provided there exists 5 > such that C{y) < C{y) for all y € Fj with y{a) = A and y{b) = B 
and \\y — y\\ < 6. 

Definition 150. A function ry € Fj is called an admissible variation provided r] ^ and 

77(a) = 77(6) = 0. 

From now on we assume that the second-order partial derivatives Luu, L^v, Luw, L^^, 
L„„, and L^^„ exist and are continuous. 
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First order optimality condition 



The next theorem gives a first order necessary condition for problem ()7.14p . i.e., an Euler- 
Lagrange type equation for the fractional /i-difference setting. 



Theorem 151 (The /i-fractional Euler-Lagrange equation for problem ()7.14p ). //y € Jt is 
a local minimum for problem, ()7.14p . then the equality 

L^im) + /.A^(,)L„[y](t) + aKL^mit) = (7.15) 

holds for all t € T*^ with operator [■] defined by [y]{s) = (s,y'^(s),aA"y(s),sAJjy(s)). 

Proof. Suppose that y{-) is a local minimum of £[•]. Let rj{-) be an arbitrarily fixed admissible 

variation and define a function $ : f — ir-mr , ii— n 

V M-)\\ ' \\m-)\\ / 

(7.16) 



He) = jC[y{-) + STii-)]. 
This function has a minimum at e = 0, so we must have $'(0) = 0, i.e., 



b ^ 



Lu[ymv^(.t) + LM{t)aKv{t) + L^ivmh^bVit) 



At = 0, 



which we may write equivalently as 

r.p(fe) 



hL.M{t)v''{t)\t=pib)+ / L^[m>f{t)At+ / L,,[y](t),A^7Kt)At 

rb 

+ / L^[y](t);,Afr?(i)Ai = 0. (7.17) 

J a 



Using Theorem 11461 and the fact that 77(a) = 77(6) = 0, we get 



b ppib) 

LMit)aA''hVit)At = I ( ,A^(b) (LM) it) ) v''it)At 



(7.18) 



for the third term in (j7.17p . Using (|7.9p it follows that 

rb 
L^[y]{t)hAlr]{t)At 

L4y]{t)ihA^"vit-J^h))^At 



Lu 



^(--l), 



/^^^^"(^ - '^h) - f(^7^(^ + ^h- amr'm 



(7.19) 



At 



£ L4y]{t)hA;^b)V''it - ^/i)At + ^0^^ j\h + vh- a{t))^r'^L^[y]{t)At. 



95 



CHAPTER 7. DISCRETE FRACTIONAL CALCULUS 



We now use Lemma 11451 to get 

f L^[y]it)hA~;i^V^it-uh)At 

J a 



L^mt) 



h'r^^{t) + 



T(v + l) 



a{t) 



{s + uh-a{t)i;;~^\^{s)A 



At 



nb 

/ /i^L^[y](i)7?^(t)At 

Ja 



+ 



ni^+'^)Ja 



Pib) 



Li,, 



Ja(t) 



At 



(7.20) 



/i^L^[^](i)r?^(t)At 



+ 



r(j^ + 1) Ja 



(u-l) 



ri'^it) / {t + uh-a{s)yr 'L^[y]{s)As 



At 



ri^{t)aAj;''{L^[y]){t + uh)At. 



We apply again the time scale integration by parts formula (|1.7|) . this time to (|7.2U|) . to 
obtain, 
rb 

la 



r,^{t)aA-^^''{L^[y\){t + uh)At 

■Pib) 



rP(o) 

/ r]^it)aA-''{L^[y]){t + i^h)At 

Ja 

+ (7?(6) - r?(p(6))),A-- {L4y]) (t + z^/i)|t=p(fe) 



[v{t)aA~''{L4y]){t + uh)] 



t=p{b) 

t=a 



Pib) 



r,%t){aAl-{L^[y]){t + vh))^At 



A^. (7.21) 



+ r?(6),A^'^ (L,^[y]) (t + z^/i)|t=p(fe) - ??(p(&))aA;^'^ (L^[y]) (t + vh%^p^,) 
r]{h)aAl'' {L^[y]) {t + i//i)|j=p(b) - 7?(a)aA-'^ (L^[y]) (t + vh)\t=a 



Pib) 



int)aAi{L^[y]){t)At. 



Since r]{a) = r]{b) = we obtain, from ()7.20p and ()7.2ip . that 

fb _^ rpib) 

/ L,.,[y](t),A;(^)r?^(t)At = -/ r?'^(t),Af (L^[y]) (t)At , 

./a ./a 

and after inserting in ()7.19p . that 

b i-p{b) 

L^[y]{t)hA^bVit)^t = / r?-(t),Af(L^[y])(t)At. 



(7.22) 



By I^TJEh and (17:^ we may write (ITTril as 

L„[y](t) + ;,A^(,) {L,[y]) (t) + ,Af (L^[y]) (t) ri%t)At = . 
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Since the values of ry'^(t) are arbitrary for t G T"", the Euler-Lagrange equation ()7.15p holds 
along y. D 

The next result is a direct consequence of Theorem 11511 

Corollary 152 (The /i-Euler-Lagrange equation - cf., e.g., |21j). Let T be the time scale liL, 
h > 0, with the forward jump operator a and the delta derivative A. Assume a, & € T, a < b. 
If y is a solution to the problem 

C[y{-)\= / L(t,y^(t),y^(i))At^min, y{a) = A, y{b) = B , 

J a 

then the equality Lu{t,y''{t),y^{t)) - {L^{t,y''{t),y'^{t)))^ = holds for all t G T''". 

Proof. Choose a = 1 and a L that does not depend on w in Theorem 11511 D 

Remark 153. If we take h = 1 in Corollary 11521 we have that 

Lu{t, r{t), Ay(i)) - AL,(t, r{t), Ay(t)) = 0, 
holds for all t G T'' (see also Theorem [3i]) . 

Natural boundary conditions 

If the initial condition y{a) = yl is not present in problem (|7.14|) (i.e., y{a) is free), besides 
the /i-fractional Euler-Lagrange equation ()7.15p the following supplementary condition must 
be fulfilled: 

- h-^LMia) + Y^^;^ ^ly + lh- af;^-'^LM{tW 

-I (t + 7/i-f^(a))i'"'^^.[y](t)At) +L^[y](a) = 0. (7.23) 

Ja{a) J 

Similarly, if y[b) = B is not present in ()7.14p (y(6) is free), the extra condition 

hLuiyMb)) + h^LMiPib)) - h'^L^mpib)) 

cb 



+ f(^ (/ (^+^^-'^w)r'^^-[^]w^* 



{p{b) + vh- a{t))f^-^^L^[y]{t)^t = (7.24) 

is added to Theorem 11511 The proofs of these facts are immediate and we don't write them 
here. We just note that the first term in (j7.24p arises from the first term of the left-hand side 
of ()7.17p . Equalities ()7.23p and ()7.24p are called natural boundary conditions. 

97 



CHAPTER 7. DISCRETE FRACTIONAL CALCULUS 



Second order optimality condition 

We now obtain a second order necessary condition for problem (j7.14p . i.e., we prove a Legendre 
optimality type condition for the fractional /i-difference setting. 

Theorem 154 (The /i- fractional Legendre necessary condition). Ify G J-x is a local minimum 
for problem. ()7.14p . then the inequality 

h^Luuim) + 2h^+'LuAym + 2h''+\iJ - l)L„^[y](t) + h^^{-f - 1)^ L„„[y]Ht)) 
+ 2/i''+T(7 - l)L,^[y]{a{t)) + Ih-'+^u - l)L,^[y]{t) + h'"'{u - lfL^4y]{t) 

+ h^''L^4y]{a{t)) + l\'L^4y]{s) (^^|i^(t + z./^ - a(s))l^~')) As (7-25) 

+ h^'L,Mit)+ I" h'L,vms)(^!f^{s + ^h-a{a{t)))i;'-'A As > 



holds for allte T^\ where [y]{t) = (t, y"(t),aAfy(t),iA^y(i)). 

Proof. By the hypothesis of the theorem, and letting <& be as in ()7.16p . we have as necessary 
optimality condition that $"(0) > for an arbitrary admissible variation i]{-)- Inequality 
<I)"(0) > is equivalent to 

/ L„„[y](i)(77'^(t))2 + 2L™[y](t)77'^(t),A^r?(t) + 2L„^[y](t)77'^(t);,Afr?(t) 

J a 

+L,Mit)iaAtvit))^ + 2L,^[y]{t)aA'^7]it)hA^^7]it) + L^^(t)UAf7?(t))2] At > 0. (7.26) 

2 

Let r € T^ be arbitrary, and choose ry : T ^ M given by 

I otherwise. 
It follows that ri{a) = rj{b) = 0, i.e., rj is an admissible variation. Using (|7.5p we get 

[Luu[y]{tW{t)f + 2L^MitW{t)aAtv{t) + L,Mit){aAt7iit)f]At 

LuuMtWit))' 

+ 2LuMitWit) [h^V^it) + ^^^^^ J\t + ^h- ais))^^-'^7^^is)As 
+ L^Mit) l^h^V^it) + ^.^-^ l\t + ^h- a(.))i^-^)^^(s)As) 



At 



■ h^Luu[y]{r) + 2h^+^LuM{r) + /i^+^L™[y](T) 
+ 1'^ L,Mit){h^V^{t) + ^,^^j\t + jh-a{s))'}J-'^rj^{s)As^ At. 
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Observe that 



h'^^Hj-l)'L,Mi^{r)) 



£^ ^..[y](i)(r(:^£(i + 7/i-^(^))i'"'^r?^(^)A.) At 



i..[y](i) /iV(i) + 



7 



air) V l-Kl + ^) Ja 



^(7-l)^A/ 



At. 



Let t G [cr2(T),p(6)] n hZ. Since 

J a 



r(7 + i) 



+ / (t + 7/i-^(s))l^"'^??^(5)As 



a(r) 



h 



7 



n{7-1) 



r(7 + i) L 

7 



r(7 + i 

7/1T 
r(7 + i) 

7/1T 
r(7 + i) 



-/i 



(t + 7/^-^(r))r^-(t + 7/i-^Kr)))J, 
■^.-1^1+4 _;^.-ir(^ + 7-1) 



(7-1) 



r(^ + i) 



r(^: 



(7.27) 



't-T 



+ 7-i)r(^+7-i) 



;¥)r(^ + 7-i) 



(¥)r(^) 

(7-i)r(^ + 7-i)' 



r(^ + i) 



27(7-1) 



(7-2) 



r(7 + 



_L^^t + jh-a{a{r)))l'-'>, 



we conclude that 



£^ L™[y](t)(^^;^-^£(t + 7/.-<7(.))t^"^)r?^(.)A.) At 



7(7-1) 



.2/^^n(7-2) 



U/™'*""rfT7Tiy*'-^^"-''<^»!' 



At. 



Note that we can write jA^r/(f) = — ^A ^i^V^it — fh) because rj{h) = 0. It is not difficult to 
see that the following equality holds: 

( 2Lu4y]{t)r]^t)hA^,r]it)At = - [ 2L„^[y](tX(t);,A-^,)r?^(t - i./i)At 

a J a 

= 2/i2+-L„^[y](r)(i/-l). 
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Moreover 

rb 



I 2L,^[y]{t)aKri{t)h^lv{t)^t 

J a 



h'r^^{t) + 



- f {s + i.h-a{t))t'^r^^is)As l At 

^) Ja{t) J J 



= 2h^+''+\u - l)L,,^[y]{r) + 2h''+''+\j - l)L,^[y](a(T)). 
Finally, we have that 

L^4y]{t){h^^^r]{t)fAt 



Lww[y]{t) 



Lww[y]it) 



i (.2^ + J^J Ja(t) 

J- (J^ + ij Ja(t) 



At 



Lww[y]{t) 



V 

u\2 , u2u+l 



{s + vh-a{t))^r^\'^{s)As 



At 



+ hL^4y]{TW - ^hy + h"'+'L^^[y]{a{T)) 



Lww[y\{t) 



V 



,(^-1) 



^vWTv) |(^ + ^^ - ^W)r ' - (^W ^y^-o 



X^ + i) 

Similarly as we did in (|7.27p . we can prove that 



w)r^}l' 



^^fl^(^+^^--(^))^'■ 



Thus, we have that inequality ()7.26p is equivalent to 



hi h^Luu[ym + 2h^+^L^Mit) + h^L^Mit) + L,,iait))ijh^ - Wf 



+ 



7(7 - 1) 



n(7-2) 



As 



la{a(t)) Vr(7+1) 

+ 2h^+^Luu,[y]{t){y - 1) + 2K<+^{v - l)L,^[y](t) 
+ 2/i^+^(7 - l)L,„^(a(t)) + /i2-L^^[y](t)(l - z.)2 + h^- L^^[y\{a{t)) 

+ 1 h'L^Jy]is)(p^^{t + uh-a{s)r-A'As\>0. (7.28) 
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Because h > 0, ()7.28p is equivalent to ()7.25p . The theorem is proved. D 

The next result is a simple corollary of Theorem 11541 

Corollary 155 (The h-Legendve necessary condition - cf. Result 1.3 of [21]). Let T be the 

time scale KL, h > 0, with the forward jump operator a and the delta derivative A. Assum,e 
a,b £T, a < b. If y is a solution to the problem 

rb 

C[y{-)]= / L(t,y'^(t),2/^(t))At^min, y{a) = A, y{b) = B , 

J a 

then the inequality 

h^Luuim) + 2hLuv[ym + L^vMt) + Lv^Mait)) > (7.29) 

holds for all t € T«^', where [y]{t) = {t,y''{t),y'^{t)). 

Proof. Choose a = 1 and a Lagrangian L that does not depend on w. Then, 7 = and the 
result follows immediately from Theorem 11541 D 

Remark 156. When /i — )■ we get a{t) = t and inequality ()7.29p coincides with Legendre's 
classical necessary optimality condition Lyy[y]{t) > (cf. Theorem [33|) . 

7.4 Examples 

In this section we present some illustrative examples. Calculations were done using Maxima 
Software. 

Example 157. Let us consider the following problem: 



1 



3 



^yi-)] = ^J (oA^yW) At^min, y(0) = 0, 2/(1) = 1 . (7.30) 

We consider ()7.30p with different values of h. Numerical results show that when h tends to zero 
the /i-fractional Euler-Lagrange extremal tends to the fractional continuous extremal: when 
h —^ (|7.3Up tends to the fractional continuous variational problem in the Riemann-Liouville 
sense studied in ^ Example 1], with solution given by 

, ^ 1 f^ dx , , 

yit) = ^ / -T ■ (7.31) 

^ Jo [(1 -x)(t -x)]4 

This is illustrated in Figure 17.11 In this example for each value of h there is a unique h- 
fractional Euler-Lagrange extremal, solution of ()7.15p . which always verifies the /i- fractional 
Legendre necessary condition (|7.25p . 
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Figure 7.1: Extremal y{t) for problem of Example 11571 with different values of /i: h = 0.50 
(•); h = 0.125 (+); h = 0.0625 (*); h = 1/30 (x). The continuous line represent function 

(I73ID- 



Example 158. Let us consider the following problem: 



>C[y(-)] 



I ioKy{t)f 



y%t) 



At 



mm . 



y(0)=0, y{l)=0. 



(7.32) 



We begin by considering problem ()7.32p with a fixed value for a and different values of h. 
The extremals y are obtained using our Euler-Lagrange equation (j7.15p . As in Example 11571 
the numerical results show that when h tends to zero the extremal of the problem tends to 
the extremal of the corresponding continuous fractional problem of the calculus of variations 
in the Riemann-Liouville sense. More precisely, when h approximates zero problem (|7.32p 
tends to the fractional continuous problem studied in O Example 2]. For a = 1 and h —?■ 
the extremal of ()7.32p is given by y{t) = 2^(1 — i), which coincides with the extremal of the 
classical problem of the calculus of variations 



£[!/(■)] 



Iy\tf 



y{t) dt 



mm . 



y(0)=0, y(l)=0. 



This is illustrated in Figure YT72\ for h = ^, i = 1,2,3,4. In this example, for each value of 



a and /i, we only have one extremal (we only have one solution to ()7.15p for each a and h). 
Our Legendre condition (j7.25p is always verified along the extremals. Figure 17.31 shows the 
extremals of problem (|7.32p for a fixed value of h (h = 1/20) and different values of a. The 
numerical results show that when a tends to one the extremal tends to the solution of the 
classical (integer order) discrete-time problem. 

Our last example shows that the /i-fractional Legendre necessary optimality condition 
can be a very useful tool. In Example 11591 we consider a problem for which the /i-fractional 
Euler-Lagrange equation gives several candidates but just a few of them verify the Legendre 
condition ()7.25p . 
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Figure 7.3: Extremal y{t) for ^T^ with 
Figure 7.2: Extremal y(t) for problem h = 0.05 and different values of a: a = 
([7:32]) with Q = 1 and different values of 0.70 (•); a = 0.75 (x); a = 0.95 (+); 
h: h = 0.5 (•); h = 0.25 (x); /i = 0.125 a = 0.99 (*). The continuous line is y(i) = 
(+); /i = 0.0625 (*). i*(l-0- 



# 


2/(1) 


5(1) 


HI) 


m) 




Legendre condition (|7.25|l 


1 


-0.5511786 


0.0515282 


0.5133134 


9.3035911 


Not verified 


2 


0.2669091 


0.4878808 


0.7151924 


2.0084203 


Verified 


3 


-2.6745703 


0.5599360 


-2.6730125 


698.4443232 


Not verified 


4 


0.5789976 


1.0701515 


0.1840377 


12.5174960 


Not verified 


5 


1.0306820 


1.8920322 


2.7429222 


-32.7189756 


Verified 


6 


0.5087946 


-0.1861431 


0.4489196 


10.6730959 


Not verified 


7 


4.0583690 


-1.0299054 


-5.0030989 


2451.7637948 


Not verified 


8 


-1.7436106 


-3.1898449 


-0.8850511 


238.6120299 


Not verified 



Table 7.1: There exist 8 Euler-Lagrange extremals for problem ()7.33p with a = 0.8, /3 = 0.5, 
h = 0.25, a = 0, b = 1, and 6 = 1, but only 2 of them satisfy the fractional Legendre condition 



Example 159. Let us consider the following problem: 

\3 



^ivi-)] 



iaA'^y{t)y + eihA^y{t)yAt^mm, y{a) = , y{b) = 1 . (7.33) 



For a = 0.8, (3 = 0.5, h = 0.25, a = 0, 6 = 1, and = 1, problem (I7:33]l has eight different 
Euler-Lagrange extremals. As we can see on Table 17.11 only two of the candidates verify the 
Legendre condition. To determine the best candidate we compare the values of the functional 
jC along the two good candidates. The extremal we are looking for is given by the candidate 
number five on Table [7?T1 

For problem ()7.33p with a = 0.3, h = 0.1, a = 0, b = 0.5, and 9 = 0, we obtain the results 
of Table 17.21 there exist sixteen Euler-Lagrange extremals but only one satisfy the fractional 
Legendre condition. The extremal we are looking for is given by the candidate number six 
on Table [7:21 
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# 


y(o.i) 


y(0.2) 


y(0.3) 


y(0.4) 


ay) 


47.25^ 


1 


-0.305570704 


-0.428093486 


0.223708338 


0.480549114 


12.25396166 


No 


2 


-0.427934654 


-0.599520948 


0.313290997 


-0.661831134 


156.2317667 


No 


3 


0.284152257 


-0.227595659 


0.318847274 


0.531827387 


8.669645848 


No 


4 


-0.277642565 


0.222381632 


0.386666793 


0.555841555 


6.993518478 


No 


5 


0.387074742 


-0.310032839 


0.434336603 


-0.482903047 


110.7912605 


No 


6 


0.259846344 


0.364035314 


0.463222456 


0.597907505 


5.104389191 


Yes 


7 


-0.375094681 


0.300437245 


0.522386246 


-0.419053781 


93.95316858 


No 


8 


0.343327771 


0.480989769 


0.61204299 


-0.280908953 


69.23497954 


No 


9 


0.297792192 


0.417196073 


-0.218013689 


0.460556635 


14.12227593 


No 


10 


0.41283304 


0.578364133 


-0.302235104 


-0.649232892 


157.8272685 


No 


11 


-0.321401682 


0.257431098 


-0.360644857 


0.400971272 


19.87468886 


No 


12 


0.330157414 


-0.264444122 


-0.459803086 


0.368850105 


24.84475504 


No 


13 


-0.459640837 


0.368155651 


-0.515763025 


-0.860276767 


224.9964788 


No 


14 


-0.359429958 


-0.50354835 


-0.640748011 


0.294083676 


34.43515839 


No 


15 


0.477760586 


-0.382668914 


-0.66536683 


-0.956478654 


263.3075289 


No 


16 


-0.541587541 


-0.758744525 


-0.965476394 


-1.246195157 


392.9592508 


No 



Table 7.2: There exist 16 Euler-Lagrange extremals for problem (|7.33|) with a = 0.3, h = 0.1, 
a = 0, b = 0.5, and ^ = 0, but only 1 (candidate #6) satisfy the fractional Legendre condition 

(IZ25D. 

7.5 State of the Art 

Discrete fractional calculus is a theory that is in its infancy. To the best of the author's 
knowledge there is no research paper providing any kind of result for the time scale T = /iZ, 
h > 0. Nevertheless, the three papers [151 UHl ES] are already published and contain some 
basic results for the time scale T = Z and also some methods to solve initial value problems 
with discrete left-fractional derivatives. None of these papers contain results involving any 
subject related with the calculus of variations. We have already two research papers \n\ [T8] 
within this topic and further investigation is being undertaken by Nuno R. O. Bastos in his 
PhD thesis under the supervision of Delfim F. M. Torres. 

The results of this chapter were presented at the ICDEA 2009, International Conference 
on Difference Equations and Applications, Lisbon, Portugal, October 19-23, 2009, in a con- 
tributed talk entitled Calculus of Variations with Discrete Fractional Derivatives. 
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Conclusions and future work 



This thesis had two major objectives: the development of the calculus of variations on a 
general time scale and the development of a discrete fractional calculus of variations theory. 
For the calculus of variations on time scales we contributed with a necessary optimality 
condition for the problem depending on more than one A-derivative (cf. Theoreml57p and also 
for the isoperimetric problem (cf. Theorem I63|) . Some Sturm-Liouville eigenvalue problems 
and their relation with isoperimetric problems were also considered (cf. Section [4. 2. ip . In a 
first step towards the development of direct methods for the calculus of variations on a general 
time scale we derived some useful integral inequalities^ and solved some variational problems 
(see Section [5.2p . As an "outsider" consequence of the herein obtained integral inequalities 
we were able to prove existence of solution to an integrodynamic equation (cf. Theorem llOSp . 
In fact, it is our purpose to study the existence of solutions to boundary value problems of 

the form 

y^\t) = f{t,y^{t),y'^{t)), t e [a,b]tf , 

y{a) = A, y{b) = B, 

and apply the possible obtained results to prove existence of solution to some classes of Euler- 
Lagrange equations (remember that the Euler-Lagrange equation is a second order dynamic 
equation) . To the best of our knowledge no such a study exists for a general time scale (we 
don't even are aware of such a study for discrete Euler-Lagrange equations). 

Remark 160. While studying integral inequalities and its applications we were able to derive 
some results which are not within the scope of the calculus of variations and optimal control. 
This work is published in the three research papers [411 1421 FIT]. 

We would also like to give answer to a deeper question: find the Euler-Lagrange equation 
for the higher-order calculus of variations problem without restriction (H) on the forward 
jump operator. We have some work in progress towards this end, namely, we obtain in |50j . 



^We remark that the Gronwall-Bellman-Bihari type inequahties presented in Section [6.41 for functions of 
two variables can be very useful in studying partial delta dynamic equations. 
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using the weak maximum principle on time scales [62], the desired Euler-Lagrange equation 
but in integral form. 

With respect to the discrete fractional calculus theory we proved some properties for the 
fractional sum and difference operators in Section 17.21 and then, in Section 17.31 we proved 
the Euler-Lagrange equation as well as Legendre's necessary condition. We also indicate, in 
Section [7.41 that as the parameter h tends to zero, the /i- fractional Euler-Lagrange extremal 
seems to converge to the continuous fractional one. This indicates that the discrete extremal 
can be used to approximate the continuous one. One of the subjects that we are willing to 
study within the discrete fractional calculus is the existence of solutions to the Euler-Lagrange 
equation (|7.15p . Note that this equation contains both the left and the right discrete fractional 
derivatives. Since we obtained so far two main results, there is plenty to do in the development 
of the theory of discrete fractional calculus of variations. For example, we want to discover 
what one gets with a general variable endpoint variational problem. One can also think in 
obtaining some criteria in order to establish a sufficient condition. Moreover, we want to find 
if it is possible to make an analogue study to that done in Chapter [5] in order to solve directly 
some discrete fractional variational problems. We think that this is a fruitful area with much 
to be done, in both of theoretical and practical directions, and therefore we would like to 
make part of the efforts to accomplish it. 

To close this manuscript it follows a list of the author publications during his PhD thesis: 

miBIlHaiiaiMlHSlSelllZlllHlllS]. 
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